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Abstract 

We consider the asymptotic behaviour of positive solutions u(t, x) of the fast diffusion equation 
ut = A (u m /m) = div (u m_1 Vit) posed for x £ R d , t > 0, with a precise value for the exponent 
m = (d — 4)/(d — 2). The space dimension is d > 3 so that m < 1, and even m = —1 for d = 3. 
This case had been left open in the general study [7] since it requires quite different functional 
analytic methods, due in particular to the absence of a spectral gap for the operator generating 
the linearized evolution. 

The linearization of this flow is interpreted here as the heat flow of the Laplace-Beltrami opera- 
tor of a suitable Riemannian Manifold (R d , g), with a metric g which is conformal to the standard 
R metric. Studying the pointwise heat kernel behaviour allows to prove suitable Gagliardo- 
Nirenberg inequalities associated to the generator. Such inequalities in turn allow to study the 
nonlinear evolution as well, and to determine its asymptotics, which is identical to the one satisfied 
by the linearization. In terms of the rescaled representation, which is a nonlinear Fokker-Planck 
equation, the convergence rate turns out to be polynomial in time. This result is in contrast with 
the known exponential decay of such representation for all other values of m. 
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1 Introduction 



In this paper we shall describe the asymptotic behaviour (as t — > oo) of a class of solutions it(i, x) > 
of the fast diffusion equation (FDE) 

(1.1) d t u = A (u m /m) = V ■ (vr-^u), m < 1, 

posecF^l for t > in the whole space, x € R d , in dimensions d > 3, and taking initial data 

(1.2) u(0,x) =u Q {x) > 0, 

where uo belongs to a class to be made precise below, in particular uq is bounded and decays at 
infinity like c |x| 2 /( 1 ~ m ) with lower order terms. Actually, since m < (d — 2)/d it is well-known that 
for initial data of the above form the weak solution exists and is unique for small times, and then 
extinguishes completely after a finite time T = T(m,d,uo), [37j- We are interested in the behaviour 
of such solutions near extinction, as t /* T. A detailed analysis of this question has been performed 
in a recent paper [7] for general m < 1 (even when m < 0), but rates of convergence could not be 



obtained for a special value of the diffusion exponent m, preciseh^j for m* = (d — 4)/(e? — 2). We 
refer to that paper for further references to the abundant literature on the topics of entropy methods, 
rescaling and rates of convergence for this type of nonlinear diffusion equations, cf. also [TU], [T2"] . 

El. US], EQ, [13, Hi, ED, M- 

The present paper is devoted to settle the asymptotic behaviour in the special case m = m*. We 
shall see that it falls out of the scope of asymptotic theory developed in the paper [7] for the rest of 
the values m < 1, both in the type of techniques and in the type of results. The clue to finding the 
stabilization rates of the rescaled orbits towards their equilibrium states in this special case relies on 

(i) Realizing that a suitable linearization of the rescaled flow can be viewed as plain heat flow in a 
suitable Riemannian manifold. This allows us to use the very detailed theory that has been developed 
for studying (the long-time behaviour of) such flows, see [3"2l 118) : 

(ii) Performing a study of nonlinear stability based on an interesting modification of the entropy 
methods of [7]. 

The paper gives precise statements and proof of these assertions. It is organized as follows: in the 
next section we shall review the needed facts about the asymptotic behaviour of our problem in the 
more general setting of variable m S R. We also introduce the family of entropies that allow to prove 
the plain stabilization result of [7], as well as the linearization method that allows to obtain rates of 
convergence when m ^ m*, when used in combination with the limit of the previous entropies. The 
failure of this approach in the special case m* is identified in as the lack of a suitable spectral gap 
in the operator analysis of the linearized problem. 

We then focus on m = m* and address such an essential difficulty. The convergence results are 
carefully stated in Section [3l We start the new work in Section [4] by a detailed analysis of the 
linearized equation, identified as a heat flow on a cigar-like Riemannian manifold. This is followed by 
the results on linearized stability. Section [5] gathers all the results needed in the comparison of linear 
and nonlinear entropies. The proof of nonlinear stability is given in Section [6l In Section [7] we revise 
the convergence for the case m ^ m* and show that our method provides a shorter proof and also a 
small improvement with respect to [7]. 



( a ) There is no restriction m > 0. The last expression represents a parabolic equation whenever u > even if m < 0. 
For m = the first expression must be replaced by Alog(u). 

( b 'ln space dimension d = 4 we have m* = 0, logarithmic diffusion. For d = 3 we deal with m„ < 0, a very singular 
case that was only briefly exposed in [7|. 
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The main difference in the asymptotic results is that convergence to a selfsimilar profile takes place 
with a rate of approach that differs in a marked way from the power rate of all the cases to =/= m*. 
The convergence is most clearly visualized below in the rescaled representation, a nonlinear Fokkcr- 
Planck equation, where it takes the form of stabilization towards equilibrium with a polynomial rate 
of approach in terms of the new time variable s. Specifically, the study is made in terms of rescaled 
variable 

(1.3) v(s,y) = (T-t)- dp u(t,x), y = ax(T-tf, s = 1 log(T/(T - t)), 

where T > is the extinction time and the constants /3, 7 and a are precisely defined in Section [jF^l 
This rescaled variable satisfies the nonlinear Fokker-Planck equation, see (|2.4[) or (|2.6p , which is better 
suited for the asymptotic analysis. The stationary profile for the latter version of the equation is given 
by the simple expression 

(1.4) V D (y) = l/(D+\y\ 2 )^/\ D > 0, 

for a suitable constant D determined by the initial data. This simple expression is handy since Vd 
and powers of it will appear as weights in some functional inequalities that are essential in our study. 
In terms of the new logarithmic time s (that goes to infinity as t — > T), the long time behaviour of the 
rescaled flow takes the form of stabilization towards the profile Vd with a power rate of convergence: 

(1.5) \\v(s,y) -V D (y)\\ L ^ {Rd) = 0(s~ 1/4 ) as s^+oo. 

This rate replaces the exponential decay formulas with respect to s of the cases to 7^ m*, that have 
been obtained in [7] . This polynomial rate in s is slower than the exponential rate in s that obtains in 
all other cases m < 1, m ^ to*. Summing up, we are in a case of what is called slow asymptotics, or 
critical slowing down, in mechanical systems and statistical mechanic, and such cases need as a rule 
special analytical methods. 

The needed assumption on the initial data is that v(0,y) be a small perturbation of Vjy(y) in a 
sense made precise by assumptions (HI') and (H2') below. Let us stress that some kind of similar 
assumption on the data is needed to obtain the asymptotic result. Actually, for data that decay 
at infinity with a slower rate than 0(|y| 2 ^ 1_m ' ) ) (i.e., with a smaller power) solutions do not even 
extinguish in finite time. On the other hand, for data that decrease with a larger power, the behaviour 
near the extinction time follows a completely different pattern that is described in the monograph 
[37]. 

We complete this introduction with some comments on related topics. Let us first recall that there 
are two other known instances of interpretation of fast diffusions as geometrical flows. The first case 
is the evolution Yamabe flow, i. e., the fast diffusion with m — (d— 2)/(d + 2), d > 3. It describes 
how a conformal Riemannian metric evolves by scalar curvature; in that case u is interpreted as the 
conformal factor of the metric raised to the power (d + 2)/4. An asymptotic study of this problem is 
made by Del Pino and Saez in [20| with exponential convergence to a separate variable solution, and 
the results are extended in [37] , In the second case we deal with Ricci flow in dimension d = 2, as 
proposed by Hamilton [25] , and then to = (logarithmic diffusion) . The asymptotic behaviour in that 
case is rather complex, cf. [16] or the monograph [37j. Both models happen in a different context, 
since they consist of interpreting the variable u in the FDE as the evolving conformal factor of a 
conformal representation, while here we consider a heat flow on a fixed manifold as the linearization 
limit of a nonlinear fast diffusion flow. They have in common the property of extinction in finite time. 

( c 'The exponent /? is essential, whereas the values of a, 7 > are just convenient. 
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Finally, we mention that a number of formulas and ideas used in the theory of Ricci flows bear a 
close similarity with developments in linear and nonlinear diffusion theory. Thus, the use of entropies 
is prominent in Perelman's study of the Ricci flow, [35] . where he introduces his functionals J- and 
W which are extensions of the Einstcin-Hilbcrt functional. He then writes the gradient flow for the 
functionals as a system of equations for the evolving metric gij and a scalar function /, which satisfies 
a backward heat equation. Strong connections exist with studies of entropies for heat equations on a 
static manifold, see for example Ni [34] and also the general references [TU [15] [29] . In a recent paper 
|33| Lu, Ni, Villani and one of the authors investigate Harnack inequalities and entropies for porous 
medium and fast diffusion equations on static manifolds that are closely related to Yau, Hamilton and 
Perelman's work, and on the other hand are close to the subject of this paper. The whole topic calls 
for further understanding. 

List of notations 

D ,Di,D*: the constants involved in Assumptions (HI), (HI'), (H2),(H2'). See Section |2~2"1 
/, /: the functions involved in Assumptions (H2), (H2'). See Section l2~2l 
^(w): the relative entropy. See Formula (12.101) . 

F: the linearized relative entropy. See Proposition 14.121 and Lemma 15.31 The argument of F can be 
both g and w (see below for the meaning of the latter quantities). 

g: the (weighted) linearization of w — 1. See Formulas (|2.14p and (12.151) . 

g Q : the metric describing the geometric interpretation of the linearized operator. See Formula (|4.3|) . 
T(w): the relative Fisher information. See Formula (|2.12|) . 

I m : the linearized Fisher information, see (|4.2[) . The index m is dropped in Section 5 for brevity. 
K(t,x,y): the heat kernel of the Laplace-Beltrami associated to g Q . See Section [4~T1 
L m : the linearized generator. See Formula (|4.ip . 

/i*: the weighted measure d/i* = y^ -m * da;. See just before Section |4~21 The L p norms in Section 4 
are taken w.r.t. 

T: the extinction time of the Barenblatt solutions and of the solutions considered. See Formula (|2.3|) 
and Assumption (HI). 

u(x,t): the solution to the fast diffusion equation. See Formula (|l.ip . 

Unit, x): the Barenblatt solutions for m > m c . Sec Formulas (|2.ip and (12.21) . 

UD,r{t,x): the pseudo-Barenblatt solutions for m < m c . See Formula (|2.3|) . 

v(y, s): the rescaled solution of the nonlinear Fokker-Planck equation. See Formula (|2.6[) . 

Vo(y): the Barenblatt profiles in rescaled variables. See Formula (|2.7p . V* (y) := Vd,(v) * s defined in 
Section O 

w: the ratio v/Vd,- See Formula (|2.9|) for the equation satisfied by w. 
Wo, W\\ lower and upper bounds for w. See Section [2~3l 

The notation || • || p denotes in principle the standard norm in L p (R d ), but starting at the end of 
Subsection UJl we will use weighted spaces and it will indicate L p (M d , d/x) with a weight /i related to 
the Barenblatt solutions. The context will always make it clear. 



4 



2 Preliminaries: rescaling, stabilization and entropy 



The fast diffusion equation with < m < 1 has attracted the attention of researchers in recent 
times, once the theory of the corresponding slow diffusion case m > 1 came to be well known. In 
the latter case the long-time behaviour of all solutions with nonnegative and L 1 data uq is given by 
a one-parameter family of explicit self-similar solutions of the form 

(2.1) U D {t,x)=t- a B D (xt- ), 

with (3 = 1/(2 + d(m - 1)), a = d/3 and profile B D = (D - k\£\ 2 ) 1 / ( - m ~ 1 '> with a free constant D > 0, 
a fixed constant k = (3(m — l)/2, and putting £ = xt^ 13 . These solutions, usually called Barenblatt 
solutions, replace the Gaussian profiles found in the long time behaviour of the classical heat equation, 
which is the case m = 1. See the precise asymptotic result in J35J Chapter 18]. 

When going over to the fast diffusion equation, the situation has been well understood in a first 
range of exponents 1 > m > m c = (d — 2)/d (the 'good' fast diffusion range); indeed, solutions of the 
above initial value problem exist and are unique, they are positive and smooth for every choice of the 
initial data in L 1 1 oc (]R d ), and even in more general cases, cf. [13]. In particular, Barenblatt solutions 
still exist, they have the same selfsimilar form though the profile looks a bit different 

(2.2) B D (0 = (D + m 2 )- 1 ^ 1 -^ 

now with k = — m)/2. This is a positive function everywhere in K d and decays at infinity 
like 0(|£|- 2/(1 - m) ), so that B D e L 1 (K d ) if m > m c . The Barenblatt solutions still represent the 
asymptotic behaviour of all solutions with nonnegative and L 1 data uq, with even better convergence 
result in relative error, cf. [8j[36]. Factors like Bjj will appear in the sequel as weights in functional 
inequalities and measure spaces. We shall use below a proper scaling to get rid of the inessential 
constant k. 

However, such a simple theory breaks down for m < m c , even if m > (which is possible if d > 3), 
due in particular to the phenomenon of extinction in finite time, cf. [3 7) . In particular, our model 
solutions cannot be continued in the same form because the similarity exponents a and (3 go to 
infinity as m goes down to m c . But for m < m c a related family of extinction solutions is found of 
the backward self-similar form 

(2.3) U D , T (t, x) = (T — t) a B D (x(T - tf), 

with (3 = l/(rf(l — m) — 2) > and a = d(3 > (just minus the formulas used before). Here, T and C 
are arbitrary positive constants and Bp is given just as in the case m c < m < 1. It is to be noted that 
Bp is no more an integrable function in K d , so we are completely away from the functional setting 
we started from. These new solutions are sometimes called pseudo-Barenblatt solutions to distinguish 
them from the original Barenblatt family. 

2.1 Rescaled flow equation 

Actually, these solutions do not possess the strong attractivity properties of their relatives for m > m c . 
In order to investigate their partial attractivity (more precisely, their rescaled stability), we have 
studied in the paper [7] the extinction behaviour of solutions with initial data close to a pseudo- 
Barenblatt solution. This is the situation in short terms: we can show that after a rescaling step we 
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obtain the nonlinear Fokker-Planck equation 

. 2 



(2.4) 



d s v = ^V„(t; m - 1 V,,v) + ■ (yv) 



in terms of the rescaled variable v(s, y) defined as 

(2.5) v(s,y) = (T-t)- d0 u(t,x), y = ax(T-tf, s = 1 log(T/(T - t)). 

Here T = T(u ) is the extinction time of the solution, (3 = (d(l — m) — 2) _1 and we will choose the 
free constants a,7 > to be a 2 = 7 = (1 — m)/3/2. Note that s(0) = and s(t) — > 00 as t — ► T. 
This means that whenever we use as T the actual extinction time of the solution u, then v is globally 
defined, for y G R d and < s < 00. With such choices equation f|2.4|) takes the convenient form 



(2.6) 



d S V = Vy ■ (V' 



V y w) 



-V? 



1 — m 

This is a convenient choice since the stationary states are now given by 
(2.7) V D (y) = (D + \y\ 2 )- 1 /^- m \ D > 0, 



which is just the profile Bo of (|2.2p without the undesired constant fc. We end this paragraph by 
noting that for m = m* the exponent in the above stationary profile is —1/(1 — m) = —{d — 2)/2, 
so that Vd(j/) decays at infinity like 0(\y\~( d ~ 2 '), i.e., like the stationary fundamental solution or 
harmonic potential. This is one of the reasons that makes to* special. 



2.2 Stabilization Result 

In paper [7] we have shown stabilization of solutions of equation (|2.6p towards one of the stationary 
profiles Vd for initial data that are not very far from Vd to start with. We can write the assumptions 
on the initial conditions in terms of either uq or vq. The assumptions on uq are 

(HI) ito is a non-negative function in L 1 1 oc (R d ) and that there exist positive constants T and Dq > D\ 
such that 

EW(0, x) < u (x) < U DuT (0, x) Vie m d . 
(H2) There exist D« G [D 1 ,D ] and f(\ ■ \) e L 1 ^) such that 

\u (x)-U D „ T (0,x)\<f(\x\) Viel d . 

In the case to < m c under consideration here, (HI) implies in particular that the extinction occurs at 
time T. Moreover, when m > to* (H2) follows from (HI) since the difference of two pseudo-Barcnblatt 
solutions is always integrablc. For to < to* this is no more true, and (H2) is an additional restriction. 

In terms of vo, conditions (HI) and (H2) can be rewritten as follows. 

(HI') vq is a non-negative function in L 1 1 oc (l£ <i ) and there exist positive constants Do > D\ such that 

v Do ( y )<My)<v Dl (y) Viei d . 

(H2') There exist D* G [£>i, A)] and /(| ■ |) G L 1 ^) such that 

\vo(y)~V D ,(y)\<f(\y\) V y ei d . 



G 



We point out that condition (HI') means a decay for large y of the form 

vo(y) = \y\- 2/{1 - m) (l-c(y)\y\- 2 ) 

with c(y) bounded above and below away from zero. Moreover, (H2') imposes a stronger decay 
condition for to < m*. Notice we can take f(\y\) = T~ d/3 f{\y\/ 'aT 13 ), so that they can be identified 
up to an elementary scaling. 

As a starting point for our asymptotic study we state the result of [7] about the convergence of v(t) 
towards a unique Barenblatt profile. 



Theorem 2.1 (Convergence to the asymptotic profile) Let d > 3, m < 1. Consider the solu- 
tion v of (|2.6[) with initial data satisfying (H1')-(H2'). 

(i) For any to > to*, there exists a unique £>* £ [Di, -Do] such that J m d(v(s) — Vd,) da; = for any 
t > 0. Moreover, for any p £ (g(m), oo], lim^oo J Rd \v(s) — Vd, \ p dy = 0. 

(ii) Form < to*, v(s) — Vd, is integrable, J„d(v(s) — Vd, )dy = 1^(^(0) — Vn,)dy andv(s) converges 
to Vd, in L p (K d ) as t — > oo, for any p £ (1, oo] . 

(iii) (Convergence in Relative Error) For any p £ (d/2, oo], 

(2-8) lim \\v(s)/V D , - 1|| = . 



For simplicity, we write w(s) instead of y i— > v(s,y) whenever we want to emphasize the dependence 
on the time s. The exponent q{m) is defined as the infimum of all positive real numbers p for which 
two Barenblatt profiles Vd 1 and Vd 2 are such that \Vd 1 — Vb 2 | belongs to L p (R d ): 

d(l - to) 
2(2 — to) 

We see that q(m) > 1 if to £ (0, to»), q(m t ) — 1, and g(m) < 1 if to > 7«». In case to > to*, the value 
of £>* can be computed at s = as a consequence of the mass balance law J Rd (vo — Vr>, ) dx = 0, and 
then the conservation result holds for all s > as is proved in the paper [7]. On the other hand, in 
the case to < to* the mass balance does not make sense, but -D* is determined by Assumption (H2'). 
In this case, the presence of a perturbation of Vd. with nonzero mass, does not affect the asymptotic 
behavior of the solution to first order. 



2.3 Relative error, entropies, and linearization 

The deeper stabilization analysis of equation (|2.6|) leads to an interesting connection with a family of 
Poincarc-Hardy functional inequalities. In this way, we obtain stabilization rates that are exponential 
in the new time s, which means that they are power-like in the original time. The exponent to* appears 
precisely as the only exponent for which the linearized analysis based on Poincare-Hardy inequalities 
fails and the corresponding rates are not obtained by that method. We shall prove below that the 
linearized analysis when m takes the special value to* leads to a different functional framework and 
the actual rates are different, and actually slower. 

In any case, the approach and the use of entropies starts in the same way. Let v be a solution to 
the rescaled Fokkcr-Plank equation (|2.6|) . and let V* = Vd, be the Barenblatt solution mentioned 
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in Theorem 12.11 We pass to the quotient w(s,y) = v{s,y)/V*{y). Notice that w — 1 = (v — V*)/V* 
is the relative error of v with respect to V*. Notice also that, by straightforward calculations, our 
running assumptions imply that Wo < w < Wi, where Wo = (D*/ Do) 1 ^ 1 "™^ < 1 and W\ = 
(ZVD x ) 1/(1 - m) > 1. 
The equation for w reads 



(2.9) 



d s w = -^V 



,,m— 1 



-v? 



In terms of w, we define the relative entropy 

1 



(2.10) 



1 - m J R d 



(w - 1) - —(w m - 1) 



V?dy 



Strictly speaking, we are assuming that a time s > is given and then we get J-(w(s)). In terms of 
v, when m is sufficiently close to one, it can be derived as E(v) — E(V*) where 



(2.11) 



E[v] =: 



vV™- 1 - —v r - 



(for m farther away from 1 both E[v] and ^[Vi] become infinite and only the expression for the 
difference, the relative entropy, makes sense). We also introduce the relative Fisher information 



(2.12) 



I[w] 



- 1 



V*wdy 



-VI: 



v dy 



(again, we should have written T[w(s)]). By differentiation in time and using the equation, we get 

dF[<u;(s)] 



(2.13) 



ds 



-I[w(s)] , 



Vs > 0. 



For a detailed proof of this time derivation, we refer to Proposition 2.6 of [7]. 

We now introduce the linearization idea in [7J that allows to treat the long-time behaviour of w. It 
consists in writing the relative error in the form 



(2.14) 



w{s,y)-l=eg(s,y)Vt m (y) 



where the choice of weight V} m is crucial. After a brief formal computation we obtain the differential 
equation for g that is implied by (|2.9[) in the limit e — > 0: 



(2.15) 



d s g = K m_2 V • (KVff) • 

implies that w — > 1 as s 



oo, the factor e will not be 



Actually, since Theorem 12.11 formula 
needed in the actual linearization step. 

Our next task is to study this linear flow; then, we shall have to relate the actual nonlinear flow to 
its linearized approximation. But let us point out that we will not need to prove the convergence of 
solutions of the original problem to solutions of the linear problem, the analysis is rather based on the 
relationship between the two linear quantities, entropy and Fisher information, associated to equation 
(|2.15[) . and the close similarity of these linear quantities and the previously defined nonlinear ones. 
These facts plus (|2.13p produce the desired convergence result. 
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In the cases m<l,m / m*, a suitable functional setting was found where the functional inequalities 
of Hardy-Poincare type corresponding to the linear flow implied the existence of a spectral gap. 
According to more or less standard theory, existence of such a gap implies exponential decay rates (in 
s) of the norms and entropy of the solutions of the linear flow. A delicate analysis of comparison of 
entropy and Fisher information between the linear and nonlinear flow allowed finally to transfer the 
result about decay rate to the original nonlinear flow. See full details in [7J. 

The problem arising when to = to* is the absence of spectral gap. We shall prove below that this is 
essential, in fact the actual rates are not exponential but power-like in s. This is related to the heat 
kernel behaviour of the operator appearing in (|2.15|) . V r " i ^ 2 V ■ (V*Vg), acting in a suitable weighted 
Hilbert space. Details will be given in Section [H where a sharp power-like decay for the heat kernel is 
proved using a most fortunate coincidence, i. e., the representation of the linear semigroup as the heat 
flow on a conformally flat Ricmannian manifold. It will also be proved that no Hardy-type inequality 
can hold for the quadratic form associated to the generator, so that it is hopeless to use the same line 
of reasoning of paper [7J . 

3 Statement of the main results for m = m* 

We are now ready to state our main results. We use the notations v(s), w(s) instead of v(s, y), w(s, y) 
and u(t) instead of u(t, x) when the dependence on time is stressed. 

We prove convergence of v(s) to the appropriate Barenblatt profile in several senses. More precisely 
we prove quantitative bounds on the convergence in suitable L p norms, on the convergence of moments, 
and on the uniform convergence of all derivatives. Convergence takes place with the same rate of the 
linearized case. 

Theorem 3.1 (Convergence with rate to the asymptotic profile) Consider a solution v of the 
equation (|2.6[) such that vq satisfies (H1')-(H2') and fix some sq > 0. Then, the entropy of the quotient 
variable satisfies 

(3.1) f[m(s)] < Ks- 1/2 V s > s • 

for some K = K(vo,Sq). As a consequence, for any j? € [0, there exists a positive constant K$ 
such that 

(3.2) ||M>( S )-V^J|| 2 <^- 1/4 Vs> So . 

The analysis of the linearized equation indicates that this rate should be optimal. We also have 
convergence without weights in suitable L p and C 3 spaces with the same rates, where we use interior 
regularity theory for parabolic equations: 

Corollary 3.2 (i) For any q € (l,oo], there exists a positive constant K(q) such that 

(3.3) \\v(s)-V D ,\\ q <K(q)s- 1 ^ V a > s ■ 
(ii) For any j 6 N there exists a positive constant Hj such that 

(3.4) \\v(s)-V D JcHw«) <^-s~ 1/4 V S >.s - 
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These power-decay results are in contrast with the exponential rates obtained in [7] for — oo < m < 1 
and m ^ m*. Rescaling back to the original space-time variables one gets the following result which 
can be called intermediate asymptotics. 

Corollary 3.3 Consider a solution u of (jl.lj) with in = m*, with initial data satisfying (H1)-(H2), 
and extinction time T. For t sufficiently close to T and for any q £ (1, oo], there exists a positive 
constant C such that: 

\\u(t) U D , (t)||, < C (T t)°^ log (T/(T - i))" 1/4 . 
with cr(oo) = d(d — 2)/4, and o~{q) = o~(oo)(q — l)/q for q < oo. 

We also obtain a quantitative bound on the decay of the relative error of v(s) with respect to Vb»- 

Corollary 3.4 (Decay of Relative Error) Consider a solution v of (|2.6p such that vo satisfies 
(HT)-(H2') and fix some sq > 0. Then for any q G {d/2, oo] and all e > there exists a positive 
constant C q such that 

(3.5) ||«(«)/Vb.-l|| 9 <C,*-T £ Vs>s a . 
If q = d/2 there is a positive constant C such that 

(3.6) \\v(s)/V D ,-l\\ d/2 <Cs-^ Vs>s . 
Finally we also have, for all j € N, that there exists a positive constant Cj 

(3.7) ||«(*)/Vb.-l|| oi(Rd) <C*- i i £ V s >.s . 

Notice that, besides having a quantitative bound, we have some other improvements on Theorem 12. II 
first because the value q = d/2 is now allowed and because convergence of C J norms is also dealt 
with. The constants involved depend also on m, d, Dq, D\, but also on the solution at time sq 
through the relative mass (conserved along the evolution) and through the uniform bound cq on the 
ratio f Rd \Vv(y)\ 2 V D (y)dy/\\v\\l < c . 

4 Analysis of the linear case 

We address now a central topic of the paper, i.e., establishing of the long-time behaviour of the 
linearized flow in the still open case with exponent m* = (d — 4)/(d— 2). The clue to our study of 
the linearized flow in this case is to interpret it as the heat flow of the Laplace-Beltrami operator 
of a suitable Riemannian manifold (M, g), with a metric g which is conformal to the standard M. d 
metric. Studying the pointwise heat kernel behaviour allows to prove Nash and log-Sobolev inequalities 
associated to the generator. Such inequalities will later on allow us to study the nonlinear evolution 
as well, and to determine its asymptotics, which will be shown to proceed with the same rate of 
convergence as the linearized one. Since the study can have independent interest, we replace g by v, 
y by x, and s by t throughout the section to conform to more standard notations. 
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4.1 Linear equation and geometry 



Given m < 1 and D > 0, wc consider the operator given on C£°(R ) (d > 3) by 



fill L V - (D + | 2 .|2^(2-m)/(l-m) v . f ^ ^ = TA 



m-2 
D 



V • (VbVu) 



We recall that for to = m* the following holds: 1/(1 - to) = (d- 2)/2, and y^ 2 (a;) = (1 + |x| 2 ) d / 2 . 
We have dropped the index * from to simplify the notation, since the particular value of D has 
no role here. We shall think of this operator as acting on the Hubert space H m = L 2 (M. d , d/z) with 
d/i = V D ~ m dx. To define it more precisely we construct the quadratic form 

(4.2) l m[v] = I {D ^ {l _ m) dx = 1 1 V,(x)l V D (x)dx, „ eC T(R-). 

Then, I m is closable in H m (for a quite general result implying the validity of the above assertion see 
e.g. [TS], Section 4.7). We denote again by — L m the unique nonnegative self-adjoint operator in H m 
associated with its closure. In fact L m has the above explicit expression (|4.ip on smooth compactly 
supported functions. There is a particular value of to for which the above operator can be seen as 
the Laplacc-Bcltrami operator of a certain Riemannian manifold (M, g), as we shall show. This in 
particular will imply (since M turns out to be complete) that L m is essentially self-adjoint on C^°(M) 
by a result of Calabi (see e.g. [18], Theorem 5.2.3). Consider indeed the following manifold, denoted by 
M, given by W 1 endowed with the Riemannian, conformally flat metric defined, in Euclidean (global) 
coordinates, by 

(4.3) g a (x) = (D + \x\ 2 )- a I, 

where I is the Euclidean metric and | • | is the Euclidean norm. We denote by /z go the Riemannian 
measure, by |g Q | = det(g Q ) the determinant of the metric tensor, by V Q the Riemannian gradient 
and by A Q the Laplace-Beltrami operator, defined on L 2 (// ga ), associated to the given metric. 

Lemma 4.1 The Laplace-Beltrami operator A a coincides with L m , precisely when a = 1 and m = 
:= [d — 4)/(d — 2), both as concerns its explicit expression (in Euclidean coordinates) and as 
concerns the Hilbert space it acts on. 

Proof. We notice that for the above choice of metric we have 



V\g^j ={d + MT Qd/2 , <?M = {D + \x\ 2 ) a d^. 

Then we have that the Dirichlct form associated to A a is given, on test functions, by 
J<*(v) ■■= J g a (V a v, V Q v)d^ gQ = J ^ y/\g a \ {x)g % i (x) ^J. ^. dx 
{D + \x\ 2 ) { - da/2)+a \V e v{x)\ 2 dx 



(4.4) 



where V e is the Euclidean gradient and the summation convention is used. Then we notice that the 
conditions that identify A a with L m : 



-(2-m)/(l-m) 



V\iM = (D + \x\ 2 ) 

v / ig^)^'(x) = (d + \x\ 2 r 1/(1 - m) s ij 
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force a, m to be related by {da/2) — a = 1/(1 — m) and da/2 = (2 — m)/(l — m). This is equivalent 
to a = 1, m = — 4)/(d — 2) = m* as claimed. □ 

We shall now compute, in the case discussed in the above Lemma, the Ricci curvature of (M,g a ). 
Hereafter we shall drop the index a, since we always choose a = 1. We put D = 1 for simplicity 
without loss of generality. 



Lemma 4.2 Then the Ricci curvature of(M,g a= i) is given, in Euclidean coordinates, by 

~ (d - 2)\x\ 2 + 2{d - 1)~ 



(4.5) 



(d - 2)xjXj 
13 ~ (1 + |x| 2 ) 2 



(i + N 2 ) 2 



Sij, 



where we write Ric = (Rij). In particular Rio on M, such lower bound cannot be improved, and 
Ric is bounded on M. Actually, Rij(x) = 0(\x\~ 2 ) as \x\ — > co in the transversal directions and it 
behaves as 0(|x| -4 ) in the radial directions. Finally, the scalar curvature is given by 



(4.6) 



R=(d-1) 



2d+ (d- 2)\x\' 



We will postpone the proof of these formulas to appendix Al not to break the flow of the exposition. 
It immediately follows that the symmetric tensor Ric is positive; indeed, given £ £ M. d , we have 



The boundedness of Ric is clear from its explicit expression. Note that for rf = 2we are dealing with 
an Einstein metric, Ric = kg (actually, it is Hamilton's cigar soliton to the Ricci flow, [551 EDi but 
for d > 3 it is not. 

Let us continue with the asymptotic analysis of the flow. By a celebrated result of Li and Yau [32] , 
the heat kernel K(s, x, y) of the Laplace-Beltrami operator of a complete Ricmannian manifold (M, g) 
with nonnegative Ricci curvature is pointwisc comparable with the quantity 

1 d 2 (x,j,) 



Vol[B(x, 

where d(- , ■) is the Ricmannian distance in (M, g), B{x,r) is the Riemannian ball centered at x and 
of radius r and Vol is the Ricmannian volume. More precisely, 

Corollary 4.3 For all small positive e there exists positive constants Ci,c 2 such that 

C-l(e) d 2 (x,y) Cnle) d 2 (x.y) 

e -T4^7)r < K(t,x,y) < e'W^F- 



Vol[B(x,y/i)\ - Vol[B{x,Vt)} 

for all x,y G M, t > 0. 

We recall that the Li- Yau bounds require completeness, a property which clearly holds for the manifold 
we are considering. We use the notation a A b = min{a, b}. 
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Corollary 4.4 The heat kernel satisfies the following properties: 

-)- 

(4.7) ^><* 

K(t,x,x) < — Vt > l,Vx G R d , 

where fi f=a /2 means that there exists two constants ci, C2 > smc/i t/iat C1/1 < / < C2/2 near to- 

t— >to 

Proof. First notice that 

rM 1 



d(0,x) = / dt 
where |x| is the Euclidean length, so that d(0,x) ~ log|x| for large |x|. Hence, 

Vol(B(0,i?))= / VW\dx= f 7JTTT1W dx 

Jb(O.R) Jd(0,x)<R [D + \x\ z ) d / 2 

d-1 

dr ~ ci?. 



R-.+00 J r<e R (D + r 2 ) d / 2 R^+oc 

Proceeding similarly, one shows that d{xo , x) w log r^K for large \x\ and hence that Vol(i?(xo, i?)) ~ 
c(i? + log I aio I ) for large i? and, say, \xo\ > 2. The short time behaviour is clearly locally Euclidean, 
with a weight depending on x given by definition by D + |x| 2 . rj 

Remark. The above corollary extends, for the present choice of the parameter m, the result of [T8] , 
Th. 4.7.5, in several respects. In fact, in the quoted Theorem the bounds on the heat kernel are from 
above and for short time only. Notice that the short time bound in the following results matches with 
the one of [18] . One may notice that, in fact, we have proved the bound 

I 

K(t,x,x)^— Vi > l,Vx G K d , 

t2 +log(l + \x\) 

although we shall make no further use of it. 

Corollary 4.5 Each solution to the linear evolution equation dtv = L rn ,v corresponding to an initial 
datum in L 1 (]R d , (D + \ x \ 2 )(™~ 2 ) / 0--™)) satisfies the bound 



(4-8) \\v{t)\\ 00 <H{t)\\v (i \\ 1 = 



ll u o||l , + 
Ci d/2 for anyO<t<l 

C2 1/2 /or anyt>l 



where Ci are positive constants. The power of t cannot be improved for such general initial data, as 
can be seen by considering the time evolution of a Dirac delta. 

Warning: Here, the symbol || • || p denotes the norm in LP(R d 7 d/z*), where d/i. = Vg _m * (x) dx, and 
we know that V^~ m * (x) = (D + \x\ 2 )~ d / 2 . This notation will be kept in the next three sections. 
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4.2 Functional Inequalities 



We recall that I mt [v] = J Rd \Vv(x)\ 2 Vd dx on smooth compactly supported functions. The domain of 
its closure will be indicated by Dom (/ m »). 

Corollary 4.6 There is a family of logarithmic Sobolev inequalities 



(4.9) / v 2 log — d M * <£/ m >] +/?(£) |MI^ 

Jr" \\\vhJ 

valid for all v £ Dom(/ m> ) n L 1 (M d , d//*) n L°°(R rf , d/z*) cmd aZZ positive e, where /3(e) = c — | loge 
/or s < 1, /3(e) = c — i loge /or e > 1, and c is a suitable positive constant. 

Proof. We have ||w(s)||oo < Cs -1 / 2 ||«o||i for large s. Interpolating between such bound and the L°° 
contractivity property (valid since I m<t is a Dirichlet form) shows that ||w(s)||oo < Cs Hx'o || 2 for 
large s. Similarly, ||u(s)||oo < Cs~ d / 4 1| || 2 for small s. The validity of such ultra-contractive bounds 
for the solution of the linear evolution considered is known to be equivalent, by [18j . Example 2.3.2, 
to the stated logarithmic Sobolev inequalities for the initial datum uq if it belongs to Dom(J mi> ) (~l 
L 1 (]R d , d//*) n L°°(M. d , d/x*). At this point the evolution has no role anymore and to avoid confusions 
we choose to write v instead of Uq in the statement. [ 

The next consequences we draw involve the recurrence of the semigroup considered. 

Corollary 4.7 The semigroup {T s } s > associated to L mr is recurrent. In particular, L rHr does not 
admit a (minimal) positive Green function and the manifold (M. d ,g a= ^i) is parabolic. 

Proof. It suffices to note that a semigroup {T s } s >o is, by definition, transient, iff J^°T s vds is a.e. 
finite for all v € L 2 (R d , d/i»). This of course does not hold in the present case because of the s -1 / 2 
behaviour for long times of the heat kernel. [ 

Corollary 4.8 There is no bounded, strictly positive, [i^-integrable function h such that 

I \v\hd^<I m ,[v] 1/2 

for all v £ Dom(/ m J. 

Proof. The existence of a function h with the stated properties is equivalent to the transience of the 
semigroup at hand, by [23J, Th. 1.5.1. □ 

Corollary 4.9 There is no bounded, strictly positive, (i*-integrable function h such that for all v € 
Dom (7 m „ ) 

/ v 2 hd[i* < I m , [v]. 

jR d 

Proof. Since h is assumed to be intcgrablc so that h d/i* is a finite measure, that we can normalize to 
1, one would have by Holder inequality that 

1/2 



/ \v\hdn,<([ v 2 hdfi,) <(/ m ,H) 1/2 

JR d \JR d / 



14 



for all v € Dom(7 m ,), contradicting the above result. □ 

Remark. The above results prove that Hardy-type inequalities relative to the Dirichlet form I mr and 
to a strictly positive integrable weight h cannot hold, even if h is required to be bounded. This shows 
that the strategy of [7J , which relied heavily on the validity of Hardy-type inequalities and allowed to 
deal with the case m ^ cannot be adapted to the present situation. 

The ultra-contractive bounds discussed above can also be related to the validity of Nash inequalities 
for J„, t . In fact we prove now some inequalities of that type in weighted Sobolev spaces which will be 
very important when dealing with the nonlinear evolution. Such inequalities play here the role that 
Hardy-type inequalities played in the case m ^ m* studied in [7J, cf. also Section [Jj The following 
crucial result is a purely functional inequality which is proved using the linear evolution only, but will 
turn out to be the key point for the study of the nonlinear evolution as well. 

Proposition 4.10 For all v such that Im,[ v \l\\ v \\\ — c o f or some c > 0, the following Gagliardo- 
Nirenberg inequality holds true: 

(4.io) hi! <Ki m .[v] 1/3 \Mt /3 , 

for all v <G L 2 (R d , d^t*) l~l Dom(/ m J, where the positive constant K depends on cq, and diverges as 
c — » +oo. 



Proof. To get the claim, first interpolate between the bound ||w(s)||oo < #( s )H u o||i an d the L 1 
contraction property to get ||w(s)j|2 < ^( s ) 1//2 |l w o||i- From this starting point we can use a known 
argument, cf. |18j . and we briefly recall it for the sake of completeness. In fact, use the semigroup 
property and the fact that I mr [v(s)] is nonincrcasing as a function of ,s to write 



H(s)\\v \\i >(«(«),«(«)) = («(2s),t)o) 



(vo,vo) 



> {v ,v ) - 2sl m ,[v }. 



dX 



Therefore, 
(4.11) 



Mis < 2sI m ,[v ]+H(s)\\v \\l 



It would then be easy to minimize the r.h.s. of the latter formula should one have H(s) = cs~ a for 
all s > 0. The fact that H(s) has such form with different powers of time when s is small and when s 
is large forces us to proceed as follows. Assuming that I mt [vo] and ||wo||i are not zero, the right hand 
side takes the value infinity both as s — > and s — * 00 hence there is a minimum for one or several 
intermediate values of s. We want to take a particular value of s that almost minimizes the above 
formula, and we want that value to correspond to the range of not small s where H(s) = C2S 
Since we assumed that I mt [vq\I\\vq\\\ < cq for some cq > 0, we consider the 1-parameter quantity 



'1/2 



Ml 

Im» N] 



2/3 



and observe that, trivially, 



s„ > 1 



a > c, 



2/:! 
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We choose a accordingly (so that it is bounded away from zero) and plug the corresponding s a into 
(|4.1ip . noticing that for s > 1 we have H(s) = C2S~ 1 / 2 ; with these choices (|4.1ip becomes 

(4.12) |kll^<Al Vo ||^ m >o] 1/3 
with K = K(a) = a + C2aT 1 l 2 . This concludes the proof. □ 

4.3 Mass conservation for the linear flow 



We introduce here the calculation of "conservation of mass" for the linear semigroup. As usual we 



put Vp~ m ' dx = d^„ 



Lemma 4.11 The following property of mass conservation holds true for every nonnegative v £ 
L x (d/x*): 



(4.13) 



-j- / v d/z* = 0. 



We give two proofs of the result, first a quite direct one and then a proof relying on the special 
geometric nature of the linear flow. 

First proof. We use the specific form of the weights involved for a direct calculation, first for an 
initial datum belonging to L 1 (d/i») n L 2 (d/i»). Choose a test function ipR supported in the Euclidean 
ball B2R with ipn = 1 on Br(0). Let t > to > and compute, for any such t: 



(4.14) 



d 

~dt 



vip R d/x„ 



Vu • VlprVd dx 



< 



'R<\x\<2R 

|Wd| \V(p R \vdx- 



-R<|x|<2ft 



\Aip R \vV D dx < 



R<\x\<2R 



c(m, d) 
R 2 



v Vd dx 



_R<|a;|<2,R 



because |VVd| < Co(m, d)Vrj / R since it is not restrictive to assume that |V<p^| < Ci/R and |A<^| < 
C2/R 2 whenever R < \x\ < 2R. By Holder inequality we obtain that 



v V D dx < ( f v 2 V 2 D m ' dx\ ( f V™> dx\ < sr R 2 

R<\x\<2R. \JR<\x\<2R. J \J R<\x\<2R J 



We obtained that 



1/2 



since we let Er := I f R< \ x \ <2R 1,2 ™* anc ^ ^ * s eas y to cnec k that J R< \ X \ <2R Vjd* ^ x — c iA 4 . 



_d 

dt 



vipR dfii 



< Ci£r 



and we notice that e R — > as R — > 00, a fact which holds because v £ L 2 ( d/i*). This proves that 



(4.15) 



lim 



v(ti)ip R dii* 



v(t )ip R d/i* 



< ci lim e R (ti - t Q ) = 

R— >oo 



for any < to < t\. We can use dominated convergence in the left-hand side, since the Markov 
property implies v(t) £ L 1 (d/i*) for all t > 0. This yields the claim for strictly positive times and for 
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initial data belonging to L 1 (d^„) n L 2 (d/x*). We can then reach t = using the strong continuity in 
L 1 (d// ! ») of the evolution, and consider general data in L 1 (d/i*) by approximation. 

Second proof. We can also use a general argument involving conservation of probability on manifolds 
with curvature bounded below. Let {T t } t >o be the semigroup associated to the Laplace-Beltrami 
operator of the manifold considered. Then {T t } t >o is a Markov semigroup and in particular it acts on 
all L p spaces (p e [1, +00]), it is contractive on any such space and it preserves positivity. We have 
shown that the Ricci curvature of M is bounded. An application of [TS], Theorem 5.2.6 then shows 
that {T t } 4 >o preserves the identity: T t l = 1. From this, conservation of the L 1 norm for data v > 
follows. In fact, with the notation v(t) — T t v and using the fact that the adjoint of T t when seen as 
acting on L 1 is T t itself but seen as acting on L°°, we have: 



\v{t)\\i = sup 

fteL~|/i|<i 



sup 

heh°°,hG[0,l] 



(T t h)v d/i : 



sup 

h6L~he[0,l] 



(T t v)hd^ 



(T t l)vdiJ, t 



Mi- □ 



4.4 Linear case. Entropy Method 

The behaviour of the heat kernel of the linear evolution considered and the L 1 contraction property 
allow to notice that for all t > to 

|| u (t)||i<Hi)||i||«Wllc»<ci!^|i. 

Notice that the above bound is sharp. In fact, consider the solution corresponding to the Dirac delta 
at ceo, namely v(t,x) = K(t, x, xq). Its L norm then satisfies, using the symmetry of the heat kernel 
and the semigroup property: 

= / K(t,x,x ) 2 dfi* = / K(t,x,x )K(t,x ,x)diJ,* 
= K(2t, x , x ) ~ ct~ 1/2 for large t 

It is easy to get the same result by entropy methods. Although this is not necessary in the present 
case due to the previous calculations, this will serve as a model for the strategy of proof used in the 
nonlinear setting, and will make already apparent the role of the Nash inequalities proved before. 

Proposition 4.12 Let F(t) = \\v(t)\\?,. Then F(t) < ctr 1 / 2 for all t > t . 

Proof. First consider nonncgative data. Having shown that the L 1 norm of such solutions is con- 
served and, moreover, using the fact that I mt [v(t)] is decreasing as a function of t, we get that 
Im, [ v {t)\l\\ v (t)\\\ < c o for all positive t and for some Co > 0. We are then allowed to use (|4.10p with 
r = 2, s = 1, so that we have 

dF(t) T . ... F 3 F 3 

—7— = -Im, [V(t)\ < -C 4 = -C- -j 

dt \\v(t)\\$ \\vq\\1 

Thus we get, integrating the above differential inequality: 

F( t )<c« 
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The same decay holds true for all L 1 data, since we may write — (i>o)- < wo < (vq)+ and use the 
order preserving property of the evolution and the decay bound already proved for nonnegative (or 
nonpositivc) solutions. In fact, denoting by v±(s) the time evolved of (fo)±, we have first that, by 
comparison, — v-(s) < v(s) < u+(s) and v 2 (s) = v+(s) + v_(s). This, together with the above decay 
property for nonnegative solutions ||w±(s)||2 < c|j(wo)±||i s_1 ^ 2 , implies that 

*■(*) = \Ht)\\l - h-mi + < £ M^h^M±Ml < ci^oii?*- 1 / 2 . □ 



5 Nonlinear Entropy Method 



Once the linear flow has been examined and its behaviour described, we prepare the way for the proof 
of convergence with rate of the nonlinear flow via a new version of the entropy-entropy dissipation 
method. We shall use the entropy and Fisher information introduced at the end of Section [2 The 
results of this section hold for any m < 1, but the main interest is in employing them for the case 
m = m* as is done in the subsequent section. From now on we revert to the notations for space, time 
and flow variables introduced in sections [1] and [2j Thus, w = w(s,y). 



5.1 Comparing linear and nonlinear quantities. The Fisher information 

We have to prove the basic inequalities that relate the linear and the nonlinear quantities of the entropy 
method. We start the analysis by a new inequality between linear and nonlinear Fisher information, 
then we recall a Lemma of [7] which compares linear and nonlinear entropy We shall write V* instead 
of Vd, ■ We put 

(5.1) X m {w] = f 

which is the (nonlinear) Fisher information. It can be linearized, as done in [7], by letting w = 
1 + egV* m and taking the limit as e — > 0. We obtain the linearized form of the Fisher information, 
that takes the expression of the Dirichlet form typical of the linearized equation 

(5.2) I m [w]= [ \V{w-l)V?- 1 \ 2 V*dy= [ |V 5 | 2 Kdy 

the relation between g and w is g = (w — 1)V"J™ _1 ; it is not restrictive to let e = 1 in the sequel. The 
next Lemma compares in a quantitative way I m and T m . This is a first attempt that will be improved 
subsequently for m = m* and m ^= m*. We drop the subindex m from both quantities for brevity. 



-VI' 



V*wdy, 



Lemma 5.1 Let < Wq < w < W\ < +oo, be a measurable function on R ; with Wq < 1 and 
W\ > 1, and assume thatX(w) < +oo. Then the following inequality holds true 

(5.3) I[w] < kil[w] + k 2 f 5 4 K 4 - 3m dy 

for any m < 1, where g = (w — l)V^ m_1 , k\ = 2VF 1 3_2m and k 2 depends only on Wi, Wq, m and d. 
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Proof. We have w — 1 = gV} m . We first re- write the Fisher information (|5.ip in the following way: 

2 



X\w] = 



(5.4) 



where we have defined 



V 



1 



m — 1 



-VT 



V*wdy 

:= / \\7(A(w)(w-l)V: n - 1 )\ 2 V^wdy 7 

JK d 



A{w) 



w 



m — l 



1 



a(w) 



(m — l)(w — 1) u> — 1 
It is easy to check that A(l) = 1, A(w) > 0, and that A(w) — > 0, when w — > oo. Moreover, 

„,,m-2 _ 



(5.5) 



A'{w) 



w — 1 



< 



since the function a(tu) = (ui m_1 — l)/(m — 1) is concave in w, so that its incremental quotient A(w) 
(taken in w = 1) is a non- increasing function of w. If we let Wo < w < W\, with < Wo < 1 and 
1 < W\ < +oo, we then have the bounds 



(5.6) 



o'(Wi) < \A(w)\ < a'(W ) = W™ 



We shall also need estimates for |A'(u>)| for Wo < w < W\ as above, and in fact it is easy to check 
that A' (I) = (m — 2)/2 and that A' is bounded away from zero. Letting now tobea function, noticing 
that w — 1 = gV}~ m and that ()5.5|) can be rewritten as (w — l)A'(w) + A(w) = w m ~ 2 , we get 



(5.7) 



V [A(w)(w- l)K m_1 ] =V[4(«))j] = A(u;)Vg + A'(u;)[Vu;]g 
= A(™)V<? + A» [V(l + gV}- m )]g 
= A{w)Vg + A'(w)gVt m [V$] + A'H.9 2 [W„ 1_m ] 

- [iH+i'W^-iJlVs + i'Hff^VK 1 -'"] 



ID 



m-2 



\7g + A'(w)[\7V}- m ]g 2 



Now we can use this equality in (|5.ip to get: 



IN 



> 



|A(w)V5 + [VA(w)]5| 2 V*wdy 
\w m - 2 Vg + A'H [VK 1_m ]ff 2 | 2 V.wdy 
\Wg\ 2 V,w 2 ^ n - 2)+1 dy- f g 4 \A'(w)\ 2 \VVi- m \ 2 V*wdy 



>lwi m ~ 3 [ \Vg\ 2 V*dy- [ g*\A'(w)\ 2 \VVi- m \ 2 V*wdy, 

where we have used the inequality |a + 6| 2 + |6| 2 > (1/2) |a| 2 valid for any a,b e E, and the bounds 
Wq <w<W\. Thus, we have 



f \VgfV*dy<-^T\ 

jR d W 1 



1 



IT" 



2m-3 



fl 4 A'H V K 



1— 771 1 
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We next remark that the weight 

4M 2 1 



VK 1_m (y) K(y) 



(D + \y\^(D+\y\2)^ 

< 1 I = t = W ^ 

{D + \y\z) 3 (D+\y\2)^ (D+\y\^ 3+ ^ 

is intcgrable whenever (d — 6)to > id — 8). Notice that when to = m* = (d — 4)/(e? — 2), £/ie weight is 
integrable. We conclude by estimating |^4'| < &o so that 

%] = / l v -9| 2 K dy < 2Wl~ 2m T[w\ + AW^'^ko / g 4 V 4 ~ 3m dy. 

JR d JR d 

This concludes the proof. □ 



5.2 Evolution properties of the Fisher information 

We now describe some further properties of the Fisher information Z[u>(s)] as a function of time, such 
as the fact that it is uniformly bounded for large s and it goes to zero as ,s — > +oo. We prove a 
new differential inequality for the Fisher information. Indeed, by Proposition 2.6 of [7J, it is easy to 
see that the Fisher information is finite almost everywhere and is the time derivative of the entropy 
almost everywhere so that: 

(5.8) F[w(s )] - F[w(s)] = [ XKO]d£ 

J s 

taking the limits s — > oo and sq — > 0, recalling that < JF[w(0)] < +oo, < JF[w(s)] — ► as s — ► +oo, 
we can conclude that Z[u>(s)] is integrable (and nonnegative) on (0, +oo). 



Proposition 5.2 In addition to the running assumptions, suppose that v(0) — Vd > 0. Then, the 
following differential inequality for the Fisher information holds true 

(5.9) dl[w(s)} ^ KiX[w(s)] _ K2l i[ w ( s )] 

as 

the constant K\ depends on m,d,So, and the constant K2 depends on m,d, the relative mass f R d{vo — 
Vn)dx, Wq and W\. Moreover, I[w(s)] goes to zero as s 00. 



Remark. The time derivative of the Fisher information is usually controlled by means of the Bakry- 
Emery method (cf. [5]) that allows to obtain spectral gap estimates. Such an estimate cannot hold 
when to = to* since there is no spectral gap. The above proposition can be viewed as a substitute 
for the Bakry-Emery method and gives a solution for asymptotic estimates in applications with no 
spectral gap. 

Proof. The proof is divided in several steps. We use the notation il = (- — m _^ p — ^ in this section 
for brevity. Note that for large s, O is uniformly bounded and |f2| < Vf^~ \v — Vd\- 
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• Expression of the derivative. We first perform a formal time-derivative of X, but in this case 
is convenient to write it in terms of v and Vd instead of w, where we recall that w = v/Vd 



dJ _ _d 
ds ds 



\Vn\ 2 vdy 



(5.10) 

= 2/ VO^V « m - 2 ^ )vdy + / |Vft| 2 ^dy = (A) + (B) 
7 R i V ds / ds 



Now wc treat the two terms separately. 
• Estimating the term (A): We have 



(5.11) (A) 



2 [ uVQ ■ V ( v m - 2 ^- ) dy = -2 ( V • [uVfi] dy. 
Je" V d V ds 

Using the equation, v s — V ■ 0, we get 

(5.12) (A) = -2 / [V • («Vfi)] 2 u™- 2 dy = -2 / [V • (wVQ)] 2 — %^ dy. 

JR d JR d il v 

Then wc have 

o [/ Rd |V.(,V»)||»|dy] 2 J-j^V^dy] 2 

^ -W _2 f n2„,2-mJ„ ^(«) 



(5.13) 

where in (i) we have used the Holder inequality 



J Rd Wv 2 -™ dy ~ w J Rd Wv 2 -™ dy 



" 2 c 2 / Kd ( V (0)-Vb)dy " ' 



hi [J/ndM] 2 
J h 2 dfi 

while in (ii) we use integration by parts, after noticing that |a||6| > ab. The point (iii) relies on 
fact that the difference between two Barenblatt solutions behaves like V D ~ m , and on the fact that 
V Do < v(t) < V Dl , so that 

dy = f ( wm ~ 1 ~ 1 ) 2 V $m- V-m d 

7 R d V m - 1 / 

„m— 1 _ i ^ 2 
777 — 1 



^ dy 

(a) < max{lf ™- 2 , W^" 2 } f \w - 1\ 2 V% dy 

JR d 

r v m 

= m & x{W™- 2 ,Wr 2 } \v-V D \ 2 -%-dy 

JR d V D 

(b) < c max{W m - 2 , Wr 2 } [ \v - V D \ 2 V™- 2 dy 



(c) < c oCl ma X {W ^ 2 , W?- 2 } f \v- V D \V 2 - m Vg- 2 dy 

JR d 

(d) =C2 [ \v- Vd| dy = c 2 f (v - V D ) dy = c 2 f (v(0) - V D )dy, 
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where in (a) we have used (|5.6p . namely 

W?- 2 \w-l\ < 



(w r ' 



(to - 1) 



< W™~ 2 \w- 1|, 



while in (b) we have used Vd > cqVd 1 and in (c) we have used \v— Vd\ < ciV^ - " 1 . In the last step (d) 
we have used hypothesis (H2') together with the fact that v(0) — Vd > and conservation of relative 
mass, proved in Proposition 2.3 of [7]. 

• Estimating the term (B). We shall use the celebrated Benilan-Crandall estimates g], that for 
solutions to the un-rescaled FDE dtu = Au m /m read 

u{t, x) 



dtu(t, x) < 



(1 - m)t 



for any t > 



if to < 1, even for to < 0. We perform the scaling to the Fokker-Plank equation, like in section [21 so 
that the Benilan-Crandall estimates read 



(5.14) 

d s v(s,y) < 



< 



[d(l-m)-2](l-m) 
2 

[d(l-m)-2](l-m) 
«i(m,d, so)«(s,y), 



+ 



1 



d(l-m)-2 (i - m ) ( e «(i-™)Ki-™)-2]/2 _ i) 
d 1 



d(l - TO) - 2 (1 - TO) (e«o(l-m)[d(l-m)-2]/2 _ ^ 



u(s,j/) 
u(s,y) 



if s > so > 0. We remark that K\ — > +oo when so but this will not be a problem. We finally 
estimate (B) 



(5.15) 



(B) = [ \\7Cl\ 2 d s vdy <C(m,d,s ) [ \\/n\ 2 vdy 

JM. d JR d 

= Ki(m,d,so)X[w(s)]. 

This calculation is formal and has to be justified, but before wc do that let us draw a first consequence. 



• Integrating the Differential Inequality. We obtained a closed differential inequality for 
the Fisher information X[w(s)] =T{s) 



dl(s) 
ds 



KlJ(s) + K 2 1 2 {s) < 



which is of Bernoulli type and can be estimated explicitly Indeed, the exact solution on (si,s) C 
[0, +oo) of the Bernoulli ordinary differential equation Z'(s) — KiZ(s) + KiZ 2 {s) = is given by 



(5.16) 



Z(s) 



a Kl(s— Si) 



< 



< C 



«1 



Z-i + £ e»itt-«)K 2 dd K 2 e"itt-0 d^ K2 



if s > si + 1 := so, for a suitable c > 1 which can be taken to be arbitrarily close to one by choosing 
si large enough. By comparison it is clear that X(s) < Z(s) < ck\/k2, provided I(so) < ^(so)- 
Therefore, for all c > 1 and < sq < s, 



(5.17) 



I[w(s)] < ^i- 

K 2 
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The constant «2 depends on m,d, the relative mass f R d{vo — Vb)dx, Wq and Wi; the constant K\ 
depends on m, d, So and ki — * +00 when sq —* 0. 

• Justification of the calculation. The differentiation of 1 performed above contains calcula- 
tions that are not justified in principle since they involve differentiations and integrations by parts in 
integrals over the whole space that are not justified a priori. Therefore, we introduce a cutoff function 
Cn{y) for the integrand of I and define 

We assume that ( n has value 1 if \y\ < n, value of \y\ > 2n, and |V£ n | < 1/n, |A£ n | < 1/n 2 . Then 
we have 

(5.18) ^ = (A^) + (B n ) 

ds 

and the two terms are as before but for the cutoff factor. There is no problem with (B n ). But (A n ) 
produces extra terms that we must control. Indeed, 

(A n ) = 2 J vwn ■ v (V"" 2 -^) £dy 

= 2 [ V • („Vfi) v™~ 2 ^£ dy 2 [ v™- 1 ^ (Vn ■ VC 2 dy. 
J R d ds J Rd ds 

When we replace dv / ds by its value according to the equation, the first of the two terms of the last 
expression becomes 

(5.19) (A nl ):=-2[ \V-(v^)\ 2 v m - 2 Cdy(=-2 [ \v s \ 2 v m ~ 2 C dy] . 

JR d \ JR d / 

which has a convenient negative sign. We now perform integration by parts on this term, an operation 
that is now perfectly justified, and we get much as before: 



The numerator of the last term is larger than | J Rd V • (i>VQ) QC« dy\, hence 



2 



(Anl) < ~«2 



/ (V-(«Vfi))n£d|, 



with the notation that we have used above. Let us calculate the integral: after integrating by parts, 
it gives a term as before plus a term where C, 2 , is differentiated, as follows: 

f v I Vft| 2 £ dy + 2 [ ufiVfl ■ C„VC„ dy = (X' n ) + (X^) 
Jm d Jm d 

The first term is (X' n ) = T n , as before, while the new term, (X' n r ) 1 can be tackled as follows. We 

1 /2 

separate by Holder a factor like X n (but we only need to integrate in the annulus R n = {n < \y\ < 2n} 
so it goes to zero as n — > 00) and we still have another factor: 

f v\n\ 2 \VCn\ 2 dy<C f vlV^-V^lVU'dy^C f ~ 2 dy 

JR d J R„ J R„ n 
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and this tends to zero as n — ► oo for m > m* . For m < m» we calculate differently, 



f vlV™- 1 ~ ^"TlVCnl 2 ^ < C / ^V|z; - Voldy < C ( \v- V D \dy 7 

that goes to zero as n — > oo, but in a uniform way we only know that is bounded a priori. In any 
case, raising to the square we get an estimate of the form 

(5.20) (A nl )< -k' 2 II + 4T, 

with constants uniform in n. 

• We now consider the other new term 



(A n2 ) = -4 / v m -\ n ^- (Vfi • VCn) dy. 



Use Holder and the numerical inequality 2ab < ea 2 + b 2 /e, to separate a term like (A n i) in formula 
(|5.19| . with a factor e which is convenient to be absorbed by the term already existing that has a 
negative sign in front, as we have remarked. We are then left with another factor of the form 

(A n22 )<c[ v m \vn\ 2 \v( n \ 2 d y ~c [ ^ivfr™- 1 -^- 1 )! 2 ^. 

JM d Jn<\y\<2n n 

Since v m ~ 1 (s,y) ~ |y| 2 as \y\ — > oo, and we get an equivalent expression 



(A n22 )<C f \Vn\ 2 v dy < CX 2n . 

Jn<\v\<2n 



I n<\y\<2n 

If we had CX„ instead of CT 2n we would have ended. The estimate for (B n ) has no problems. 

• To solve the difficulty we take a little detour. We integrate the inequality obtained so far for dX ra / ds 
to get the integrated inequality: 

(5.21) In{s2)-In{si)<k / X n ds + C l 2n ds 

J si J si 

But the right-hand side is bounded above by the integral 

/•S2 

(5.22) (C + k) Ids 

J Sl 

and this is known to be bounded by the relative entropy. Moreover, since the integral X ds is 
finite, for every e > there exists a s £ such that 

/oo 
I ds < £. 

We conclude that I n {s 2 ) — X n (si) < e when s E < s\ < s 2 . Combining this half continuity with the 
integrability of X„(s) < 1{s) given by (|5.23p . we obtain by an easy calculus lemma that 

(5.24) I„(s) < Cie 
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for all s > 2s £ , with C\ uniform in n. We conclude that I(s) = lim n X„(s) is bounded for all large 
times and goes to zero as s — > oo. 

• Coming back to the differential inequality satisfied by T n we have proved that I' n < c\I — c^T^ . 
Integrating this differential inequality in time between s\ and S2 with si < s 2 sufficiently large we get 
In{s2) — In(si) < c\ J s ^ I{ s ) ds ~ c 2 J s * In{s) 2 ds so that, passing to the limit as n — > +00 and using 
both monotone convergence and the boundedness of / as a function of time we get /(S2) — I(si) < 
Ci J s ^ I(s) ds — C2 J Si 2 I(s) 2 ds, which is an equivalent form of our statement. Q 



5.3 Comparing linear and nonlinear entropies 

The quantitative comparison of linear and nonlinear entropies concludes the preliminary results needed 
for the nonlinear entropy method. Under Assumptions (H1")-(H2"), the relative entropy is well 
defined. 

Lemma 5.3 (An equivalence result) Let m < 1. If w satisfies (H1")-(H2"), then 
(5.25) J>L<T [W] <J>L. 
We recall that F[w] = J Rd \w - 1\ 2 V% dx. 

The short proof of this result has been given first in [7] but we repeat it here for reader's convenience. 
Proof. For a > 0, let (j> a (w) := j^— [(w - 1) - (w m - l)/m] — a (w — l) 2 . We compute <j)' a {w) = 
— ^ [1 - w" 1 ' 1 ] -2o(k;- 1) and tj%(w) = w m ~ 2 - 2a, and note that cf> a (l) = (f>' a (l) = 0. With 
a = W[ n ~ 2 /2, <f)'a is positive on (Wo,Wi), which proves the lower bound after multiplying by Vg l 
and integrating over R d . With a = W m ~ 2 /2, </>" is negative on (Wq,Wi) which proves the upper 
bound. [~J 

Equivalently, we may write 

FH F\w] ^ n F\w] F\w] 
5.26 Sr*— < < F\w\ < 77— k < 

V / r»TT7-2— m — n .... L., 2-m — l J — o:„r I „.J2— m, — 



2W(~ m ~ 2sup|u>| 2 - m ~ 2val\w\ I - m ~ 2W<f 



5.4 The entropy bounds a suitable L p -norm 
Lemma 5.4 Let m < 1. If w satisfies (H1")-(H2"), then 
(5.27) |k-l||'SSr (Ed) <D m F[w], 

where D m is given at the end of the proof. 



Proof. We first state some inequalities between Barcnblatt solutions with different constants. Con- 
sider 

^ = -T^—[D + \y\T^ = ~t~~ — ^d" 1 < . 

oD 1 — m 1 — to 

Hence, for any < D\ < Dq 

Do - Pi 2 _ m D - J>i 2 _ m 
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Moreover, it is easy to see that if < D\ < Dq 

v * m (y) = s 5TTBF " ^ " lv) 4 + " + v -™ <9) ' 

The above inequalities prove that |tu— l\ m ' ( 1_m ) is bounded by a multiple of V™ . Indeed, by hypothesis 
(HI') we have that < D < D* < D±, so that Vd —V*< v(s) — V* < Vd x — V*, and since w = v/V* 



\w-l\ 

Thus, 



v(s) - V* 



< \V Dl -Vbol < J>o-£>i Vg/ 

K 1 771 K 



Ik-Ill 2 :^ . = / \w-lftw-l\**dv 



1 - 777 y v ^1 

Remarks, (i) The estimate proves that w — 1 £ L 2+ 1-™ (R d ), whenever the initial entropy is finite, 
since we know, joining inequalities (|5.25| and (|5.27j) : 

(5.28) \\w(s) - l||'tS (R(i) < < 2^ m ^ 1 2 -™^[ U ;( S )] < 2A„W?-* n .F[uto], 
since the nonlinear entropy is decreasing in time. Moreover, we have also proved that 

(5.29) |Ks) - lll^j^ < 2D m Wt m nw(s)} 
and we shall show below that entropy goes to zero as s — > +oo . 

(ii) As an easy consequence, letting w — 1 = (v — V,)/V, and using the fact that V* < C, we obtain 

(5.30) \\v-V4 2 j £% 0ii) <D m F[w]. 

(iii) For m = m* we have 2 + j^— = d/2. 

5.5 Comparing linear and nonlinear Fisher information 

With the above remarks we can improve on Lemma HTT1 that compares the linear and nonlinear Fisher 
information: 



Proposition 5.5 Under the same assumptions of Lemma \5.1\ we have 
(5.31) I[g] = / |V.g| 2 V* dy < kiX[w] + kzF 1+a [w] 

for anym<l, where g = (w - 1)KT~\ o = 2/[d + 2 + 777/(1 - m)) > 0, k\ = 2W?~ 2m , and fc 3 > 
is given at the end of the proof. 
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Proof. We estimate the second term of the inequality of Lemma 15.11 in the following way 
f g ±vt 3m dy= [ (\w-l\V m - l ) A Vt 3m dy= [ (\w - l\)\ m dy 

JWl d JR d JR d 

<!k'-l|| 2 oo / (\w-l\) 2 Vrdy = \\w-l\\ 2 00 F[w] 



Now we recall the interpolation inequality (|7.7[) with j = 

(5.32) ||/||l-(r*) < C d \\f\\%? m \\f\\f~ p 

then we apply it to / = w — 1 and we let p = 2 + m/ (1 — m). We get: 

i 

d / 9 | m \ d + p 

(5 33) \\w-l\\ L oo m < C p4 \\w - l\\ d c V m {\\w(s)-l\\^^L^ d) J 

< C p , d Nh {2D m W^ m T[w])^^ := fc 3 ^N CT / 2 

where a — 2/[d + 2 + m/(l — m)] > for any m < 1 and we used inequality (|5.28[) and the fact that 
\\w - lll^Ld) < A^i by Theorem O Thus we have proved that 



/ 9 4 Vt 3m dy<\\w- F[w] < k^[w]' 



The expression of k% is then 

where \\w - l||^ RtJ) < M u D m = D °Z°* D {2° 1 (l + 37) , and C P:d is the constant of the 
interpolation inequality 17.71 with j = and p = 2 + m/(l — to). Q 

Remarks. (i) The above proposition holds for any to < 1 and allow to conclude that I(s) — > 
as s — > +00, since we already know that both T(s) (cf. Proposition 15. 2p and J-(s) tend to zero as 
t — > +00. 

(ii) When to = to*, we obtain that a = 4/ (3c?). But in this critical case we shall need another finer 
comparison between the linear and the nonlinear Fisher information that hold only when to = m„. 
namely we would like to have that there exists sq > and a constant > 0, such that 

(5.34) I[g(s)} < kO[w(a)]. 

for any s > so, where g = (w — 1)V™*~ . Unfortunately the above inequality is not guaranteed for 
all times s > sq. In the next section we will prove a weaker version of this statement, sufficient to our 
scopes, namely we will show that the above estimate (|5.34[) holds on a family of intervals [s\ k, S2,fc] 
that is sufficiently dense as s — > 00. The technical details will be postponed to Appendix A4. 

6 Proofs of the main results in the critical case 

In this section we shall always take to = to* , and we shall show that the nonlinear flow converges with 
the same rate as the linear case, cf. Section l4~4l We shall use the relationship between the entropy 
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functional T and the Fisher information Z, namely dT j ds = —X. In view of the absence of any 
spectral gap (or Hardy-Poincare inequality) inequality valid instead in the case to ^ to* , we have to 
proceed differently The Gagliardo-Nirenbcrg inequality that in the linear case give the correct decay 
of the linearized entropy in the L 2 (V„ 2_m da;)-norm, turn out to work as well in the nonlinear case as 
the previous proposition started to show. 

Proof of Theorem 13. li Notice first that ||<?||oo is finite and bounded as a function of time. In fact, 
by hypothesis (HI') we know that 



\g(s,y)\ = \w(s,y)-l\V™- 



v-V D , 



Vd, 



V, 



m — 1 



<co\V Dl (y)-V Do (y)\V™ 



m-2 



t r2 — in, t rm—2 
< C 1 V D, V D, 



Cl 



for all y € R d and all s > 0, where a arc a positive constant depending only on to, Do, D±, -D*. The 
inequality \Vd 1 — Vb | < cV D ~ m can be proved easily using the explicit expression of the pseudo- 
Barenblatt solutions (see the proof of Lemma [57^ 



Next we prove that I is bounded as a function of time. Indeed by Lemma 15. II we observe that 
(6.1) I[g] < krl[w] +k 2 f g 4 V^ 3m dy < kil[w] + k^WgW^ 

JR d 

where we have noticed that, for to = to*, V^ _3m = (l + |a;| 2 ) 2 is integrable. It has been proved 

in Proposition 15 . 21 that X is also bounded. 

By conservation of relative mass, cf. Proposition 2.3 of [7], we know that ||g(s)||i = ||g(0)||i, where 
we have used the fact that v — Vd, above is taken also nonnegative with J (vo — Vd, ) dy = M > in 
this part of the proof. This implies that the ratio I /M = I m , [<7(s)]/||<7(s)||i is bounded as a function 
of time. 

We shall use now the Gagliardo-Nirenberg inequalities of Proposition 14.101 taking v = g(t), putting 
F = 11.911^2(^2-™ dx ), I the linear Dirichlct form, and M = \\fj\\L 1 (V' 2 - m dx) : 



(6.2) 



F 3 < K\1M 



The validity of such inequalities depends on the boundedness of the ratio I m , [g(s))/\\g(s)W^ i^ v2 - m dx y 
which is ensured along the evolution, as above mentioned. 

We now prove an entropy - entropy production inequality. We obtain a differential inequality for the 
entropy by comparing it with the Fisher information X via Gagliardo-Nirenberg inequalities, 



(6.3) 



V m - 2 / | W -l| 2 ^d y 



1, 



F 4 



m-2 



KiIM 4 = K 2 IM A 



where (a) follows from (|5.25p of Lemma 15.31 while in (b) we used the Gagliardo-Nirenberg inequality 
above. 



(i) In order to continue the argument, we assume for the moment that the initial datum satisfies 
vo > Vd, and is radially symmetric so that 90 = (wo — 1)V^~ i s nonnegative. This extra assumption 
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will be removed afterwards. Under it we will prove in Appendix A4 that there is an infinite sequence 
of intervals of times [si^ , S2,k] C [2k, 2k + 2] (hence, S2,k < Si,fc+i) such that 

(6.4) I[g{s)} < ka[g{s)} for all s e |J [si, k , s 2 , fe ] 

for a constant fci that does not change along the evolution. We shall prove moreover that the length 
of each of such intervals is at least 1/2 for all k > kg, which in particular implies that 

n n \ k 

(6-5) ( S2 > fc ~ - 2 = — 2 _ ^- CS2 '™ 

k— ko k—ko 

whenever n > no is large, for a suitable c > 0. Then, recalling that X = — AT ' j ds, and using (|6 . 3[) we 
conclude that 

^ 3 < A' 2 fc 4 M 4 1 <k 5 l= ~k 5 — , 

di 

and an integration oven the interval [si,fc, S2,fc] gives 

£ G^b* " ^f) - £ S (S2 ' fc " si ' fc) - i S2 ' n • 



This implies 

1 



since the intermediate terms are such that 

1 



f(si,k) 2 T(s 2 ,k-i) 



< 



because F(s) is non-increasing and S2,fc-i < si^. The monotonicity of the function T(s) allows then 
to conclude that for all s £ [s2,fc, s 2,fe+i] 

, . 111c 

.F(s) 2 T(s 2 , kY F\s\,\Y h 

Using the fact that s < S2,k+i < s 2 ,k + 4 we get 
(6.7) T{s) < — - < 



[^"(Sl.l)- 2 +€^8] * (C + C 1S )2 



for large times s and some positive constants Co, c\. We have thus proved that the nonlinear entropy 
decays with the same rate as the linear one, when the initial relative mass is nonzero. 

(ii) Proof without extra restrictions. The arguments used above and in Appendix A4 are 
valid changing h into —h and g into — g under the same a priori bounds. Hence, the conclusion is 
valid for negative and radial initial difference vq — Vd, < 0. 

To deal with the general case where vo — Vd, is not radial or does not have a sign, we use the 
maximum principle, after writing |wo(^c) — Vd,{x)\ < /(|x|). By comparison we have v\ < v < v 2 , 
where v\ and v 2 are the solutions corresponding to initial data Vd, — / and Vd, + f resp. For the 
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corresponding w = v/Vd,, h = w — 1 and g = h(Djj, + y 2 ) a similar comparison holds. Thus, 
wi < w < W2, where w\ < 1 and W2 > 1 are the solutions with radial initial data 1 ± (//Vb, ), hence 
functions of r = \y\ and s. Same idea applies to g. Take now into account the form of the entropy 

(6.8) F[w] := I ¥(w)Vff dy, with $(w) = (w - 1) - — (w m - 1). 

We note that ^(w) is convex and has a zero minimum at w = 1. Since we have just proved that the 
decay result holds for both gi and 32 > the statement also holds for g, even if we do not assume that 
v a ~ Vb, is nonnegative or radial. fj 

Proof of Corollary 13.21 We recall the following facts proved in [7], Lemma 6.2 under the running 
assumptions, (HI) and (H2). First we have that for any 1? G [0, fES> there exists positive constants 
K$ , K2 such that 

\\\x\^(v-V D ,)\\ 2 <K^(^[w]) 1/2 . 

Moreover 

\\v- V D J 2 <K 2 (T[w]) 1/2 . 
We now recall the result of Lemma 3.6 of [7] 

(6.9) \\v(s)- Vb.||o-(H-) < H\W*)- VbJU Vt>t . 

for a suitable a € (0, 1), and we combine it with the interpolation inequality (|7.6p . with A = —ad < 

= (J, < 1/2 = v, C = 0,1/2 
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IK*) - Vb. Hoc < c IK*) - vb, IK*) - v^ ||^ < IK*) - Vb, lit IK*) - Vb, || 

where $ = 1/(2 + ad). This implies 

IK*) - Vb, Hoc < C 1 ^ 1 -") W<V(i-*) || l)(s) _ ^ || 2 < {T[w]) 1/2 V t > t . 
From Holder's inequality, 

- Vb, || g < |K«) - Vb, nr 2)/9 IK*) - Vb. II 2 / 9 < rc q (^N) 1/2 

1 /2 

for all q G [2, 00], we deduce that \\v (s) — Vb, ||q decays with the same rate as (.F[it;]) . 
If g G (1, 2), we know from Lemma 6.2. of [7] that there exists a positive constant K (q) such that 

\\v-V D , \\ q <K(q) LFM) 1/2 , 

This and the known decay of J- proves (ii). 
To prove (hi), use first (|7.7|) with the choice p = 00, i.e. 

(6.10) ||/lb(R«) < C jld ||/||^ (R(!) H/Ilt 1 

for any j G N, and the decay of the L°° norm, namely ||u — Vb,||oo < A's -1 / 4 to get 

IK*) - Vb, llojpad) < H, jS ~^+T) V s > s , 
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where in fact Hj depends on s itself and tends to zero as s — » +00, so that the bound can be improved. 
Indeed we iterate the procedure putting such bound for the C^ +1 norm into (|6.10[) to get a new bound 
for the C- 7 norm. In fact what we get after h steps is, for any fixed s > sq, j £ N: 



\\v(s) - Vb,||cw(R<«) 



< 



where the value of kh will be determined later. Notice in first place that the numerator of the above 
expression remain finite as h — > 00, for any fixed s > sq, j £ N. As for kh, by construction it satisfies 
the recursion relation 

1 3 , 1 



4(J + 1)' -"^-i + l"* ' 4(i + l) 
subtracting 1 /4 to both sides of the latter equation gives 

kh+i 



4 j + 1 V 4 

which immediately gives fc/j = | — (j+i) 4(j+i) ' ^ us P rovm § that kh — ► 1/4 as ft. — > +00. [] 

Proof of Corollary 13.41 We have proved that ^*[iy(s)] < cos -1 ' 2 and by Lemma \5. 31 we also know 
that F[w] < ciJ-"[w]. By Lemma 15.41 with to = to*, we have 

(6.11) \\w(s) - lfj* = \\w(s) < D m ,F[w(s)] < c 2 T[w{s)] < c^ 1 ' 2 

Moreover, (|5.33| yields 



(6-12) \\ W ( S ) l|| Loo(Ri) < CA T[w{s)Y^^ = c,T[w{ S )] 2 /^ < c^ 1 '^ 

Interpolating between these bounds shows that, for q £ [d/2, +00]: 

\\w(s) - l\\ q < \\w - l\\^™i\\w(s) - ifj^ < ces-itiH). 



To improve such bound we insert it in the interpolation inequality (|7.7[) and use Theorem l7.2l as well 
to get 

C 

\\w(s) - lllc-j^) < qll n fc _, 

for any g £ [d/2, +00], k > j £ N. As a function of q the exponent of s is nonincreasing, so we choose 
q = d/2 to get 

C 

IhO) - l||ci(R<«) < ■ 

S d(2 + fc) 

To optimize in k we should take k = 00, which is not allowed, so that for any fixed e > we take k 
large enough so that 

C 

IK«) ~ 1||cW(r*) < - i=r 
s a 

as claimed. Putting this bound back into (|7.7j) with j ' = and using what is known so far for the 
decay of the L p norm we get a decay of the form \\w — l\\oo < Cs~ a , a being given (with an inessential 
renaming of the free parameter e) by 



a(p, k) = 



l-e\ d 1/1 1\ pk 



d + pk 3 \ d p J d + pk 
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Maximizing a w.r.t. to p when e is small enough yields again p = d/2, so that after some calculation 
we get the exponent a(d/2,k) = ^ — d M^ fc \ ■ This proves the claim for the L°° norm and hence also 
for all h q norms with q € (d/2, +00) by interpolation. [J 



7 Proofs for fast diffusion with m m* revisited 

The previous method allows for shorter proofs of the convergence when to ^ , and at the same 
time some minor improvements of paper [7]. We recall that, in the case m 7^ m*, the spectral gap 
inequality 

(7.1) F[g(s)] <X m ] d I[g( S )} 

holds true, and the best constant is known for m < m«, since [H [7]. We also recall the result of 
Proposition 15.51 

%] = / |V 5 | 2 V D , dy < k!l[w] + k 3 T 1+ °{w] 
where, in particular, k\ = 2W\~ 2m . From these bounds we get 

T(w) <<<*> V™- 2 / \w - 1| 2 ^ dy = \w™~ 2 F 



( 7 - 2 ) < (b) \w^X m ) d I[g] 



< (c) \wr 2 K) d [kiAM + hf 1+a W] 

where in (a) we compared the linear and nonlinear entropies via inequality (|5.25[) . in (6) the above 
spectral gap inequality, and in (c) the above mentioned Proposition 15.51 We may rewrite the latter 
formula as a differential inequality: 

W 2m ~ 3 

T' + W W\ Am,d^ - k 3 T 1+a < 

so that by comparison with its explicit solution, we get 

p -fe 4 (s-So) 

F[w(s)] < — r < fc 5 e- feiS 

provided .F[u;(so)] is small enough, a property which holds for sq large enough. This gives an ex- 
ponential decay of the entropy, with a rate k^ = Wq _ ™W r 1 2 " l_3 A mi£ 2. Note that \ m ,d is the optimal 
constant in the Hardy- Poincare inequality (|7.1[) . known for m < to* since [7]. Then we can proceed 
as in [7j to show that the optimal rate is given by X m ,d- Indeed, one can substitute Wq and Wi with 
inf K <i \w\ and sup R d \w\ respectively, allowing them to depend on time. Then prove that they both 
tend to 1 when s — > +00, so that fej — > X m .d] this can be done in view of the uniform convergence of 
the relative error ||u>(i) — IIIloo^-) — > as s — > +00 together with a Gronwall-type argument. For 
more details we refer to Section 6.3 of [7]. 

Remarks 
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(i) When m = to* the above steps do not hold since we do not have a spectral gap for the linearized 
generator. This is one of the reasons which forced us to use Gagliardo-Nirenberg inequalities which, 
instead, compare the Fisher information with a power of the entropy. 

(ii) This method simplifies and complements some proofs of [7J when m 7^ m*, but also gives a more 
detailed proof of the case m < that was only briefly treated in [Jj . We finally emphasize the analysis 
of the present paper covers the case to = 0, that is logarithmic diffusion, even in dimension d = 4 
since in that case m* = so that no spectral gap holds. 

(hi) The interpolations made in the proof of Theorem 13. 21 are valid also in the case to 7^ to* and allow 
to improve the convergence rate of the derivatives, proving that the rate is always given by \ m ,d ■ We 
state here this improved version of the main asymptotic Theorem of [7J 

Theorem 7.1 (Convergence with rate, to 7^ to*) Under the assumptions of Theorem 12.11 i/to 7^ 
m* , there exists t > such that the following properties hold: 

(i) For any q £ [g*,oo], there exists a positive constant C q such that 

\\v(s)-V D J q <C q e- x ^ s Vs>s . 

(ii) For any $ € [0, (2 — to)/(1 — to)], there exists a positive constant K$ such that 

|| \x\*(v(s) - VdJ|| 2 < c- x ^ s V s > < • 
(hi) For any j G N, there exists a positive constant Hj such that 

Hs)-V D ,\\ C i m < H 3 e- X ™« s Vs > s ■ 

The constants C q , and Hj depend on sq, to, d, Vq, Dq, D\, and q, ?9 and j; sq also depends on 
Do and D\. It is remarkable that the decay rate of the nonlinear problem is given exactly by \ m ,d- 
Rescaling back to the original equation, we obtain results in terms of intermediate asymptotics, cf. 
Corollary 13.31 or Corollary 1.3 of [7J . 

Appendices 

Al. Calculation of curvatures. Proof of Lemma 14.21 

We can use well-known formulas for the Ricci tensor as a function of the metric data: 

Rij = g krn Rikj m , Rihjm = 2 ( d lj9irn + df m g k j - dl m gij ~ d^g km ) + 9n P {^kj T im + r km T ij)- 

but in the case of conformal transformation there is a worked out relation between the Ricci tensors of 
two metrics g and g in terms of the conformal factor relating them. Precisely, if g = (1/V 2 )g, where 
<p is a scalar, the formula reads as follows [5]: 

R-R = \[{d- 2)(pV 2 <p + (<pA s <p - (d - l)g(Vp, Vtp) ■ g] . 
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where R = (Rij) is the Ricci tensor of g, R = Rij is the Ricci tensor of g, V denotes the gradient, V 2 
the Hessian and A g the Laplace-Beltrami operator with respect to g. Specializing the formula to the 
case g = 6{j, so that Rij = 0, we get in coordinates 

Rij = ^2 [(<* - 2 )^> + (<P*<P - ( d - 1)|V^| 2 )%] • 
Put now gij = (1 + \x\ 2 )- 1 5 lj so that ip = (1 + \x\ 2 )^ 2 . Then d l (p = x,{\ + \x\ 2 )- 1 / 2 , 

ay = - XiXj (i + |*| 2 r 3/2 + (i + \x\T 1/2 6io, = -\ x \ 2 (i + \x\ 2 )- 3 / 2 + d(i + M 2 )- 1 / 2 . 

Applying the last formula we get after some calculations 



(d — 2)xiX i 

1T+W) 2 



(d- 2)\x\ 2 + 2(d- 1) 

(i + M 2 ) 2 



There is a clear form of these expressions when we take the particular point x = (X, 0, -,0) which 
implies no loss of geometrical generality since the metric is conformal and radial, hence invariant under 
rotations in the space. We get 

/ s ~ ^ 2(d-l) ~ ,^ (d-2)X 2 + 2(d- 1) w 

(7.3) fl 11 (S) = _L__^ ; R u (x)= [ - 1— ^1 > Vi = 2,...,d, 

and Rij{x) = for all t 7^ j. Both eigenvalues tend to zero as \x\ — > +00 with different rates. It 
immediately follows that the symmetric tensor Ric is positive; indeed, given £ € M. d , we have 



and the same is true for all x £ M. d by invariance under rotations. If one wants to visualize the 
behaviour of this manifold, it is convenient to look at Ricci curvatures given by 

~ i?(ei, ei ) 2(d-l) _ i?(e 8 ,e ? ) 2(d - 1) + (d - 2)X 2 
(7.4) n = — = — ttttt. r, - 



S(ei,ei) (1 + X 2 )' 1 5(e i>ei ) (1 + X 2 ) 

Note that the transversal curvatures tend to (d — 2) as \x\ — > 00 while the curvature in the radial 
direction behaves like O(|o;|~ 2 ). This clearly shows the difference in the behaviour of the curvatures 
in radial and transversal directions which is typical of a cigar manifold. 

Finally, the value of the scalar curvature follows from the formula R = g^Rij. Since we are in a 
conformal situation, it can be deduced in a direct way from the Yamabe formula |39j 

~ _ 4(d-l) Aw _ (d _ 2)/4 

d-2 w (d+2)/(d-2) ' witnw-s 

where g is the conformal factor, here (1 + |x| 2 ) -1 , cf. the formulas e.g. in [37], pages 211-212. In 
order to obtain the results stated in Lemma 14.21 we only need to eliminate the tildes from Rij and 

it □ 
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A2. Explicit representation of the cigar 

We give here a simple parametric representation for the cigar-like manifold (M, g). We want to 
represent it as a hypersurfacc in The radial symmetry of the metric suggests to represent such 

imbedded manifold as z = f(\y\) with variables (y, z) £ where y £ M d and z £ R, having a 

unique chart x £ M. d i— > (y, z) given by the formulas 

r = \y\ = $(g) and z = ^(g), 

where q = \x\ > 0. We fix <I>(0) = ^(0) = 0. We want the Euclidean metric in R d+1 to induce the 
metric on the hypersurface. We know that the infinitesimal length element in the radial direction 
satisfies 

ds 2 = dr 2 + dz 2 = [*' 2 ((?) + $' 2 (£>)] dg 2 



dg 2 



1 + g 2 

which implies the relation *f?' 2 (g) + $ /2 (f?) = 1/(1 + g 2 )- On the other hand, the length calculation for 
the transversal part gives 



g 2 ' 



Solving the above two equations gives 



*(*) = 9 ne) - g(2+g2)1/2 

(0) (l + g 2 ) 1 / 2 ' {e> (1 + g 2 ) 3 / 2 ' 

We see that r = <&(g) goes from to 1 as < g < oo. Analyzing the behaviour of ^(g), one concludes 
that 

^(g) as g 2 when g as 0, ^{g) ~ log Q when g^S> 1. 

This is the representation of a cigar. The point out that the transversal radius at infinity is constant; 
actually, $(f?) — > 1 as g — > +oo. 



A3. Some Technicalities 

We recall here some technical facts that we used in the proofs. First we recall Theorem 2.4 of [7], 

Theorem 7.2 (Uniform C k regularity) Let m < 1 andw £ L^ c ((0,T)xR d ) be a solution of (|2~9|) . 
Then for any k £ N, for any sq £ (0, T) , 

(7.5) sup ||w(s)|| c . fc ( K d) < +oo . 

s>s 

We also needed an interpolation Lemma due to Gagliardo [24], cf. also Nirenberg, [301 P- 126]. 

Lemma 7.3 Let X, fi and v be such that — oo < A < /i < v < oo. Then there exists a positive constant 
C\,^ yV independent of f such that 

(7.6) wmj* < c x ^\\mj^ \\mj v x v / e c(« d ) , 

where || • 1^/^ stands for the following quantities: (i) Lf a > 0, i/ien ||/||i/<r = |/| 1 ^' 7 da;)' 7 . (mJ 
J/ a < 0, /ei fc &e £/ie integer part of {—o~d) and a = \a\d — k be the fractional (positive) part of a. 
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Using the standard multi-index notation, where \rj\ = r/i + . . . + % is the length of the multi-index 
rj = (rji, . . . rjd) G we define 



l/lll/<T 



max \d'f = max sup : = \f\\ Ca(Rd) if a > 

\n\=k ,a \ v \=k x yl£Rd \x -y\ a 

max sup \d ,1 f{z)\ := H/HcW) if a = 
\m= k zef 



.As a special case, we observe that \\f\\-d/j = \\f\\c^(w L ) ■ 

(Hi) By convention, we note ||/||i/ = sup^ eK d \f(z)\ = \\f\\c°(&) = ||/||oo- 

Remark. The following special case of the above interpolation inequality (|7.6[) has been used in the 
paper: let k > j e N and A = — k/d < /x = — j/d <v=l/p. Inequality (|7.6p becomes 

(7.7) H/llcW) < 11/11?^ 

for any k > j E N and p > 0. 



A4. Complete proof of the estimates of Section [6] 

In the proof of Theorem 13.11 in Section [5] we have assumed that for every solution under the stated 
conditions there is an infinite sequence of intervals of good times [s^j. , s 2i fc] C [2k, 2k + 2] with 
s 2,k < s i,fe+i for all k, such that 

(7.8) I[g{s)} < k 4 l[g{s)} for all s e (J [s hk , s 2 , k }- 

fceN 

Recall that in view of hypothesis (H2) and the discussion made in Section [6l we may assume that g 
is radially symmetric and positive. We will also prove that the length, l k = S2.fc ~ si.fc of the intervals 
in our construction is at least 1/2 for all k > k . 

The proof of these facts is long, and would have broken the flow of the proof of Theorcm l3.ll this is 
the reason why we put it here. The main point in getting (|7.8[) consists in obtaining stronger estimates 
of the remainder term in the inequality of Lemma 15.11 than the ones obtained in Proposition 15.51 We 
restate here Lemma 5.1 for convenience of the reader: 

Let < Wq < w < W\ < +oo, be a measurable function on M. d , with Wo < 1 and W\ > 1, and 
assume that I(w) < +oo. Then for any m < 1 the following inequality holds true 

(7.9) I[w(s)} < kaWs)} + R[w(s)}, with R[w(s)] = k 2 f g(s, y) 4 ^) 4 " 3 " 1 dy , 

JR d 

where g = (w — l)V7 n_1 ; k\ and k 2 are positive constants. 

We need to control the remainder term i?[w(s)] to proceed with the asymptotic estimate. Note that 
y-4-3m j g i n tegrable for m = m*: in fact, for such a value of m we have 



3G 



Put now N(s) = N[g(s)] = |j<?(-, s)||oo, the supremum of g for fixed time s > 0. Wc know that N(s) 
is uniformly bounded in time. Then we have 



(7.11) 



R[w(s)} < k 3 [N(s)Y 



We want to estimate the decay of N(s) in time in terms of the linearized Fisher information I[w]. 
Suppose for a moment that we can prove that the remainder term is small relatively to I[w], more 
precisely that R[w(s)] < ^I[w(s)] for all large s. In that case we conclude that /[u'(s)] < 2kiX[w(s)], 
and the desired estimate (|7.8p easily follows. Hence, we need to prove that 



(7.12) 



N(s) 4 < 



K 



with K > ks, say K = 2k% . This is a most convenient estimate on the values of g. 

Unfortunately even under such assumption it is not clear that the last inequality holds at all times, 
or even at all large times. Therefore, we shall be cautious and call good times those times at which 
(|7.12|) holds with K > 2k^. The frequency and density of the intervals of such times is important, as 
the end of proof of Theorem 13.11 shows . 

We will now proceed with the proof of the existence of the time intervals stated at the beginning of 
this section. They will consist only of so-called good times. The proof is split into two parts, namely 

(i) Controlling the remainder term away from the origin. This is the part where we use the fact the 
l^o — ^* I < / f° r a radially symmetric /; 

(ii) Transforming the outer control into a control on a small ball, namely we will control s\vp r>R g{r, s) 
for small R > 0. Due to the peculiarities of the parabolic Harnack inequality, we shall prove that such 
a control only takes place for a large set of so-called good times, this sufficing for our goals. 

Part (i). The control of a radial g far away from the origin 

In the calculation that follows we drop the s-dependence for convenience since time does not enter in 
the argument. Let g(r) be a nonnegative continuous function such that g(r) — > as r — > oo and let 



Mi 



9(r) 



(l + r 2)d/2 



r d 1 dr < oo , 



and 



If 



(]_ + r 2)(d-2)/2 



" dr <C oo . 



We put the powers in a way such that it is clear that for r > 1 we are dealing with the radial case and 
we are merely asking the mass and the linearized Fisher information to be finite. Now pick a > 0, 
Ri > R > 1 and calculate 



,1+a 



(R)-g 1+a (R 1 ) = -c 



Ri 



< a 



g'{r)g a (r) dr < a 
\g\r)\ 2 rdr 



\g' {r)\g(r) a dr < a 

H 

Rl 9(r) 2a ' 



Rl i a(r) a 

\g'( r )\ r h^_dr 



R 



V - 



Now, if we assume a > 1/2, the last integral can be bounded as follows: 



11 \9(r)f 



dr < sup | <7 C 7 ™) I " 

r>R 



\9(r)\ 



dr, 



so that for a > 1/2 we have obtained 



(R)-g i+a (R 1 )<a 



sup \g(r) 

r>l 



\9(r)\ 



dr 



\g'(r)\ 2 rdr 
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Letting Ri — ► oo, and assuming that g(R) — > as i? — > oo, we get 



3 1+a (tf) < " 



sup |fif(r 



|3'(r)| 2 rdr 



Taking the supremum over R > 1 on the l.h.s. and simplifying we get: 



supg(R) 


4 

8 

< as 




4 
3 




R>1 




.A r 







|. 9 '(r)| 2 rdr 



This is a very good estimate because it says that the supremum of g A outside the unit ball is not only 
proportional to / as expected in so-called better times, but even more: it is proportional to a higher 
power of I. Now, recall that I[w(s)} — > as s — > oo. If the same could be done near r = the proof 
that every large s is a good time would be complete. 

The previous calculation can be done in the complement of the ball of radius R as small as we like 
and then g(R) will depend also on an inverse power of R, because of the presence of the factors 1 + r 2 
in the denominators of the last quantities. We now get in the last line for < R < 1 

C 



sup \g{r)\ 


4 

8 

< a 3 


\r |5W| d / 


4 
3 


\ \g'(r)\ 2 rdr 


4 
3 

< 


r>R 




Jr r 




Jr 





R s(d~i)/3 J l R 2 R- 

The estimate blows up at R = 0. Therefore, we cannot let R — > to get an estimate for N(s) for any 
x G R d . 

Justifying that g goes to zero at infinity. To conclude part (i) of the proof, it remains to prove that 
g(R) — ► as R, — ► oo . Choose R n such that 



and 



1 

- dr < 



An 2 



and define R n = minjr > R n ;g 2 (r) < -^z} ■ Indeed the set in the r.h.s. is not empty since g/r is 
intcgrable at infinity: this is not compatible with g being everywhere larger than a positive constant 
for all r > R n - Notice that R n > R n and g 2 (R n ) < l/(2n 2 ). Hence, for all R> R n : 



g 2 (R)=g 2 (R n ) + 2 / g g> dr < 



1 
2n~ 2 



9\g'\&r 



< 



1 

2^2 



„' |2 



dr 



9 dr 



< 



1 



Therefore, < g(R) < 1/n for all R > R n - The proof of part (i) is now complete. 



Part (ii). Transforming the outer control into a control on a small ball. 

In part (i) we have estimated the supremum of a radial g 4 outside a ball of radius R > in terms of 
I[g(s)] 4 / 3 , so that the problem is to estimate the supremum inside a ball as well, hopefully in terms 
of I[g(s)) 1+a , at least in the form el[g(s)]. We arc unable to prove that for all (sufficiently large) 
times. To circumvent such a difficulty we have to make use of a rather complicated argument that 
takes into account the possibility that such estimate does not hold because of possible bad behaviour 
of g at points near the origin. We begin by carefully labeling the times. We say that a time s £ [0, oo) 
belongs to the class of good times Qk, if 



N(s) 4 



sup(.g(s,y)) < 

v 



K 
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We are not claiming that some half- line [T, oo) C Q2k 3 , which would finish the proof in the simple 
way. Finally, we say that a time is very good, s G Vc, if 



(7.13) 



iV( s )<C/[ 3 ( S )] 4 / 3 



for some C > 0, in the spirit of the radial estimate away from the origin. Note that since I(s) -^Owc 
have the inclusion of very good times with constant C into the good times with any constant K > 
if s is large enough. 

Harnack inequality. The study of points near the space origin is based on classical regularity 
theory for linear or quasilinear parabolic equations in divergence form. We are going to use the 
version of the celebrated paper by Aronson-Serrin [T] . We consider the equation satisfied by the error 
function 



1 



v(s,y) 



VdM 



1. 



(7.14) h{s,y)=w(s,y) 
It can be written in the standard form 

(7.15) d s h = V-A(y,h,Vh) + B(y,h,Vh) 



In fact, starting with the equation satisfied by w, we have 
(7.16) 



d s h = d s w = — V • 



1, 



so that we can identify 



A = (h + l)V 



(h + 1) 



m — 1 



i 



i 



-V r ' 



-VI 



(h + iy 



- 1. 



-V T ' 



B = ( t + l)^.V 



(ft + 1)' 



m — 1 



-VI 



and we have to check that the structure conditions are satisfied by A and B in a compact ball with 
constants that do not depend on s. In fact, the structure conditions are satisfied in the homogeneous 
form of pQ, which means that, in the notation of that paper, the terms /, ft, = in the structure 
condition for A, B vanish. Checking this is a straightforward calculation involving also the known 
bounds on w namely < Wo < w = h + 1 < W\ . We note in passing that since we already know that 
w — * 1 uniformly in K d as s — > oo, the lower and upper bounds Wq, W\ can be taken closer and closer 
to 1 if we restrict the time to s > sq and sq is large enough. 

In any case, we conclude that the Harnack inequality has the standard form, as stated below. This 
also implies a similar Harnack inequality for g if we work on a bounded space domain, say, in £>i(0). 
We state it next. Take T > large and consider the parabolic cylinders 

Q=[r-2,T]xB 1 (0), Q 1/2 = [T-l/2,T]xB 1 (0), Q = [T - 2,T - 1] x B^O) . 
The parabolic Harnack inequality on the disjoint cylinders Q1/2 and Q is then 
(7.17) inf g(s, y) > c sup g(s, y) = c N 

Ql/2 Q 

for some positive constant c < 1 depending only on structural constants. Note that since g is contin- 
uous on Q all the above suprema are attained at some points. 
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Evolution of the maximum of h. The equation for h can be written as 



d.h = V 



+ V™~ 2 (h + l) m - l VV* • V/i 



(ui)v(M^r- 1 ) 

+ K m_3 |VK| 2 (/i + l)[(h + i)" 1 - 1 - 1] = 

= v/i ■ v f (/t + 1)m 1 1 ~ 1 K"- 1> ) + (/i + i)v ■ + rp-Vr^v/i] 



TO — 1 

+ mV: n - 2 (h + l)™" 1 W* • Vft 

+ (/i + i)p + l)" 1 " 1 - i]V • (KT" 2 VK) + K m_3 |VK| 2 (ft + i)P + i) m_1 - i] 

In particular, using the fact h is small we get 

(7.18) d s h — second and first order terms in h + C(r, h) h 

where C(r, h) < k for a suitable k independent of r and depending only on the known a priori bounds 
for h. Then, as a consequence of the Maximum Principle, cf. e.g. [T], the maximum N(s) obeys the 
growth rate 

(7.19) N'{s)<kN(s). 

In fact the function H(s,y) — N(s )e kl - s ^ So ' > is an explicit supersolution in the whole space for s > sq. 
The structure of good times. Alternative. Now we are going to prove that 

Lemma 7.4 For every time interval (T — 2,T) with T large enough there is at least a subinterval of 
length 1/2 consisting of very good times. 

Proof. We can use the cylinders Q1/2 and Q as in the previous paragraph on the Harnack inequality. 
The idea is to consider separately the two possibilities: (a) either the maximum in x of h at every 
time of the lower cylinder is taken outside the ball i?2(0), or (b) the maximum at one of such times, 
say so, is taken inside. 

In case (b), T must be a good time if it is large enough as we show next. We take s £ (T — 2, T — 1) 
as above and let Ni be the corresponding maximum in the ball. Of course Ni < N. We now apply 
the growth rate of previous paragraph to obtain 

(7.20) N(s 2 ) < CTVj < CN, C = e 21 . 

for every S2 in the upper cylinder: T — (1/2) < S2 < T. On the other hand, for every y £ Bi(0) we 
have for such S2 the lower estimate h(s2, y) > cN. We conclude that the maximum and the minimum 
at all those times are related by a constant. This is also true for the function g up to a small change 
in the constant, hence 

(7.21) 9(s2,y'u)<Cig(s2,y) 

where now y' M is the point of maximum of g in the ball -Bi(O). Now, we know that on the boundary 
\y\ = 1 there is a good estimate for g, more precisely, for such y of unit norm g(s2,y) satisfies the 
estimates that defines the very good time, and it does with a fixed constant C. We conclude that 
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g(s2,y), \y\ > 1, also satisfies such an estimate with a possible worse constant C = C\C. Since the 
estimate was true for \y\ > 1 we are done. 

Therefore, whenever T is not a good time in Qk with K = 2k%, the first part of the alternative, (a,), 
must be true. But in that case for every time in the lower cylinder we know that the maximum of h 
is taken in the exterior of the ball B\. If we look for the expression of g = h (1 + \ y\ 2 ) this also means 
that the maximum of g(s, •) at times in (T — 2, T — 1) is taken at an exterior point (maybe different). 
So all these times are very good. Recalling what was said before, they are good times in Qk if T is 
large enough. □ 

Choice of intervals of good times. We can apply the previous results letting T = 2k + 2, for 
k > fco and fc sufficiently large. The above lemma implies that there exists a subinterval [si^ , s 2 .fe] C 
[2k, 2k + 2] of length at least 1/2 made of times in Qx with K = 2k 3 . 

8 Concluding remarks and open problems 

[8J.1. The special situation has been studied for the critical exponent to* = (d— 4)/(d— 2) in dimensions 
d > 3 where our considerations make sense. Algebraic extensions have been shown to be fruitful or 
intriguing in some dynamical studies. Here, for d = 2 we formally have to* = ±oo, which is an extreme 
situation for porous medium that has appeared in the literature (for instance, in connection with the 
mesa problem), |22l [51 [3"T], while for d = 1 we formally get to* = 3, a value inside the porous medium 
range where nothing special has been shown to happen. 

[8J.2. We pose the following questions: 

• Are the rates obtained in this paper optimal for a certain class of data, as the linearized analysis 
suggests? 

• Can we prove convergence, maybe with worse rates or without rates, for more general initial data? 
we recall that for to > m c all nonnegative initial data in L 1 (R d ) are attracted towards a Barenblatt 
solution, with no rate in that generality. 

• Find an explicit optimal dependence of the constant in the asymptotic formula with respect to the 
data. 

• Assuming that we get an optimal rate of convergence, can we find a profile for the next level of 
approximation? 

[5J.3. One may wonder if the techniques used in [6l [7] for the case to ^ m*, which use Hardy- 
Poincare inequalities, work also in the case to*. We have partially given a negative answer to this 
question in Corollary 14.91 in which we have shown that no inequality of Hardy type can hold for 
the linearized Fisher information I[w\. However, one may wonder if modified versions of the Hardy- 
Poincare inequalities, with logarithmic terms added in the spirit of the classical Hardy inequality in 
R 2 , allow to solve the problem. Thus, there is a family of valid Hardy inequalities (see below) in which 
the Dirichlet form involved has a logarithmic correction, but we are not able to prove the asymptotic 
results by means of such inequalities. It is then a further open problem to see whether this path may 
lead to the goal or not. 

Proposition 8.1 Let d > 3. We have 



(8.1) 
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for any g G r D(M. d ) and for any < a < | — 1, where 

, a-l 



(8.2) 



d/i a (y) = (l + |y| 2 ) 5 [l + log(l + M 2 )f dy, 
d^(y) = (1 + M 2 ) 1 ^ [1 + log(l + |y| 2 )f +1 dy 



and 

(8.3) H a ,d = 



2{d-2) 



: (d — 2 — 2a) min {2a, (d- 2 -2a)} 
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